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ABSTRACT 


This  report  describes,  through  additions  and 
modifications  to  the  Level  15.1  NASTRAN  User's 
and  Theoretical  Manuals,  a  thermostructural  capability 
for  NASTRAN.  In  addition  to  this  new  rigid  format, 
a  set  of  two-dimensional  and  three-dimensional  iso¬ 
parametric  finite  elements  was  added  to  NASTRAN's 
finite  element  library.  The  thermostructural 
capability  consists  of  computing  a  temperature  history 
of  the  structure,  taking  account  of  such  thermal 
conditions  as  radiation,  convection,  flux,  and  heat 
generation,  and  then  performing  a  series  of  structural 
static  analyses,  using  as  part  of  the  static  loading  the 
equivalent  loads  due  to  the  temperature  distribution  at 
times  selected  by  the  user. 

This  version  of  NASTRAN  is  available  for  the  UNIVAC 
1108  and  CDC  6000  computers. 


ADMINISTRATIVE  INFORMATION 


This  project  was  funded  under  NASA-Defense  Purchase 
Request  H-84825A  issued  by  the  George  C.  Marshall  Space  Flight  Center, 
Huntsville,  Alabama.  This  work  was  performed  under  work  unit  1844-011. 


INTRODUCTION 


This  report  describes  the  thermostructural  capability 
which  has  been  added  to  the  NASTRAN  structural  analysis  program. 
The  report  consists  of  modifications  and  additions  to  the  Level  15.1.1 
NASTRAN  User's  and  Theoretical  Manuals.  The  page  numbers 
on  the  following  pages  correspond  to  the  page  numbering  in  those 
manuals.  A  brief  description  and  general  flow  of  the  new 
capability  are  given  in  Chapter  19  of  the  Theoretical  Manual. 

This  version  of  NASTRAN  is  available  for  the  UNTVAC  1108 
and  CDC  6000  computers. 


USER’S  MANUAL  ADDITIONS 


STRUCTURAL  MODELING 


The  material  property  definition  cards  are  used  to  define  the 
properties  for  each  of  the  materials  used  in  the  structural  model. 

The  MAT1  card  is  used  to  define  the  properties  for  isotropic  materials. 

The  MATS1  card  specifies  table  references  for  isotropic  material 
properties  that  are  stress  dependent.  The  TABLES1  card  defines  a 
tabular  stress-strain  function  for  use  in  Piecewise  Linear  Analysis. 

The  MATT1  card  specifies  table  references  for  isotropic  material 
properties  that  are  temperature  dependent.  The  TABLEM1 ,  TABLEM2 , 

TABLEM3,  and  TABLEM4  cards  define  four  different  types  of  tabular 
functions  for  use  in  generating  temperature  dependent  material  properties. 

The  MAT2  card  is  used  to  define  the  properties  for  anisotropic 
materials.  The  MATT2  card  specifies  table  references  for  anisotropic 
material  properties  that  are  temperature  dependent.  This  card  may 
reference  any  of  the  TABLEM1,  TABLEM2,  TABLEM3,  or  TABLEM4  cards. 

The  MAT3  card  is  used  to  define  the  properties  for  orthotropic 
materials  used  in  the  modeling  of  axisymmetric  shells.  This  card  may 
only  be  referenced  by  CTRIARG,  CTRAPRG  and  PTORDRG  cards.  The  MATT3 
card  specifies  table  references  for  use  in  generating  temperature 
dependent  properties  for  this  type  of  material. 

The  MAT7  card  is  used  to  define  the  properties  for  isotropic 
materials  for  isoparametric  elements  and  for  BAR  elements  in  thermal 
problems.  The  MATT7  card  specifies  table  references  for  isotropic 
material  properties  that  are  temperature  dependent.  This  card  may 
reference  any  of  the  TABLEM1 ,  TABLEM2 ,  TABLEM3 ,  and  TABLEM4  cards. 

The  MAT8  card  is  used  to  define  the  properties  for  anisotropic 
materials  for  two-dimensional  isoparametric  elements.  The  MATT8  card 
specifies  table  references  for  anisotropic  material  properties  that 
are  temperature  dependent.  This  card  may  reference  any  of  the  TABLEM1, 
TABLEM2 ,  TABLEM3  and  TABLEM4  cards.  ~ 
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The  MAT9  card  is  used  to  define  the  properties  for  anisotropic 
materials  for  solid  isoparametric  elements.  The  MATT9  card  specifies 
table  references  for  anisotropic  material  properties  that  are  temperature 
dependent.  This  card  may  reference  any  of  the  TABLEM1,  TABLEM2, 

TABLEM3,  and  TABLEM4  cards. 

If  MATT 7,  MATT8,  or  MATT9  cards  are  used  in  Rigid  Format  14, 

TRANGE  cards  must  be  specified  for  each  referenced  table.  TRANGE 
cards  will  cause  the  tabular  functions  specified  on  TABLEM1,  TABLEM2, 
TABLEM3,  and  TABLEM4  cards  to  be  changed  into  step  functions  to  be  used 
in  the  thermal  transient  analysis. 

The  GENEL  card  is  used  to  define  general  elements  whose  properties 
are  defined  in  terms  of  deflection  influence  coefficients  and  which  can 
be  connected  between  any  number  of  grid  points.  One  of  the  important 
uses  of  the  general  element  is  the  representation  of  part  of  a  structure 
by  means  of  experimentally  measured  data.  No  output  data  is  prepared 
for  the  general  element  is  the  representation  of  part  of  a  structure  by 
means  of  experimentally  measured  data.  No  output  data  is  prepared  for 
the  general  element.  Detail  information  on  the  general  element  is 
given  in  Section  5.7  of  the  Theoretical  Manual. 

Dummy  elements  are  provided  in  order  to  allow  the  user  to  investigate 
new  structural  elements  with  a  minimum  expenditure  of  time  and  money. 

A  dummy  element  is  defined  with  a  CDUMi  (i  =  index  of  element  type, 

1  <  i  >  9)  card  and  its  properties  are  defined  with  the  PDUMi  card. 

The  ADUMi  card  is  used  to  define  the  items  on  the  connection  and  property 
cards.  Detail  instructions  for  coding  dummy  element  routines  are  given 
in  Section  6.8.13  of  the  Programmer's  Manual. 

1.3.2  Bar  Element 

The  bar  element  is  defined  with  a  CBAR  card  and  its  properties 
(constant  over  the  length)  are  definedwith  PBAR  card.  The  bar  element 
includes  extension,  torsion,  bending  in  two  perpendicular  planes  and 
the  associated  shears.  The  shear  center  is  assumed  to  coincide  with  the 
elastic  axis.  Any  five  of  the  six  forces  at  either  end  of  the  element 
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may  be  set  equal  to  zero  by  using  the  pin  flags  on  the  CBAR  card. 

The  integers  1  to  6  represent  the  axial  force,  shearing  force  in  Plane 
1,  shearing  force  in  Plane  2,  axial  torque,  moment  in  Plane  2  and 
moment  in  Plane  1  respectively.  The  structural  and  nonstructural 
mass  of  the  bar  are  lumped  at  the  ends  of  the  element,  unless  coupled 
mass  is  requested  with  the  PARAM  card  COUPMASS  (see  Section  3.5.1). 
Theoretical  aspects  of  the  bar  element  are  treated  in  Section  5.2 
of  the  Theoretical  Manual. 
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STRUCTURAL  MODELING 


1.3.10  Isoparametric  Elements 

NASTRAN  includes  six  isoparametric  elements,  two  planar  elements, 
two  solid  elements,  and  two  surface  elements.  A  one  dimensional  element 
is  also  available  for  use  on  problems  using  isoparametric  elements.  The 
two  planar  elements  are  of  the  membrane  type  only.  The  elements  are 
defined  by  connection  cards  as  follows: 

1.  CIS2D4  -  planar  linear  element  with  four  grid  points. 

(quadrilateral) 

2.  CIS2D8  -  planar  quadratic  element  with  eight  grid  points. 

(quadriparabolic) 

3.  CIS3D8  -  solid  linear  element  with  eight  grid  points. 

4.  CIS3D20-  solid  quadratic  element  with  twenty  grid  points. 

5.  SURF1  -  one-dimensional  surface  element. 

6.  SURF4  -  linear  surface  element. 

7.  SURF8  -  quadratic  surface  element. 

Theoretical  aspects  of  these  elements  are  treated  in  Section  5.13 
of  the  Theoretical  Manual. 

All  the  isoparametric  elements  except  the  three  surface  elements 
may  be  used  as  either  structural  or  thermal  elements.  The  user  must 
include  an  APPISO  card  to  indicate  which  element  to  use.  As  thermal 
elements,  however,  they  may  be  used  only  in  Rigid  Formats  1  and  14.  The 
three  surface  elements  are  dummy  elements  which  are  used  to  conveniently 
apply  thermal  boundary  surface  conditions  in  Rigid  Format  14.  Also,  SURF4 
and  SURF8  can  be  used  to  apply  a  uniform  static  pressure  load  on  a  surface 
of  an  isoparametric  element. 

The  properties  of  the  IS2D4  and  IS2D8  elements  are  given  on  PIS2D4 
and  PIS2D8  cards,  respectively.  The  IS3D8  and  IS3D20  elements  have  no 
corresponding  property  cards.  Material  properties  for  the  planar  elements 
are  specified  on  MAT7  and  MAT8  cards.  These  properties  may  be  made 
temperature  dependent  by  including  MATT7  or  MATT8  cards.  Material 
properties  for  the  solid  elements  are  specified  on  MAT7  or  MAT9  cards 
and  may  be  made  temperature  dependent  by  including  MATT7  or  MATT9  cards. 
Examination  of  these  cards  will  show  that  the  isoparametric  elements 
may  have  isotropic  or  anisotropic  materials  and  that  the  thermal  as 
well  as  structural  material  properties  may  be  temperature  dependent. 
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Either  lumped  or  consistent  mass  matrices  will  be  computed  for  the 
structural  elements  and  lumped  or  consistent  capacitance  matrices  will 
be  computed  for  the  thermal  elements.  Lumped  matrices  will  always  be 
computed  unless  a  PARAM  card  is  presented  with  parameter  COUPMASS 
specified  (see  Section  3.5.1).  Consistent  matrices  are  usually  mandatory. 

The  element  coordinate  systems  for  the  planar  elements  are  shown 
in  Figures  13  and  14.  Symbols  Gl,  G2.  ...  refer  to  the  required  order 
of  the  connected  grid  points  on  the  connection  cards  defining  the  ele¬ 
ments.  The  angle  0  for  the  planar  elements  is  the  orientation  angle  for 
anisotropic  materials.  (For  solid  elements,  material  coordinate  systems 
must  be  set  up  when  necessary.  For  further  discussion,  see  Section  5.13.2.2 
of  the  Theoretical  Manual.) 

Element  stresses  at  grid  points  are  computed  when  requested.  The 
following  real  membrane  stresses  are  output  for  the  solid  elements  on 
request: 

1.  normal  stresses  in  the  x-,  y-,  and  z-directions . 

2.  shear  stress  on  the  x-face  in  the  y-direction. 

3.  shear  stress  on  the  y-face  in  the  z-direction. 

4.  shear  stress  on  the  z-face  in  the  x~direction. 

The  following  real  membrane  stresses  are  output  for  the  planar  elements 
on  request: 

1.  normal  stresses  in  the  x-  and  y-directions . 

2.  shear  stress  on  the  x-face  in  the  y-direction. 

Presently  no  complex  stresses  are  available. 

Also  available  are  the  average  stresses  at  the  grid  points,  as  well 
as  the  maximum  element  grid  point  stresses.  These  represent  the  normal 
output  from  Rigid  Format  14.  In  other  rigid  formats,  a  DMAP  alter 
will  be  required.  For  further  discussion,  see  Section  3.15  of  this  manual. 
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Figure  13  -  Planar  Quadrilateral  Linear 
Isoparametric  Element  Coordinate  System 


1.3-23  (7/1/72  NSRDC) 


STRUCTURAL  MODELING 
APPLIED  LOADS 

1.5.4  Thermal  Transient  Loads 

NASTRAN  includes  four  types  of  thermal  transient  loads.  For  a  full 
discussion  of  these  loads  see  Section  19.1.1  of  the  Theoretical  Manual. 

The  RADIAT1  card  defines  radiation  heat  exchange  between  a  surface  of 
an  isoparametric  element  and  its  surroundings  in  the  form 
Q(t)  =  HA(T4  -  T4) 

where 

Q  is  the  net  rate  of  radiation  heat  exchange 
H  is  the  radiation  factor 
A  is  the  area  of  the  surface 
T  is  the  temperature  of  the  surface 
T^  is  the  ambient  temperature. 

H  may  be  temperature  dependent  and  is  given  on  the  TABLEH1  card.  T^  may 
be  time  dependent  and  is  given  on  the  TABLET1  card. 

The  RADIAT2  card  defines  radiation  heat  exchange  between  two  isopara¬ 
metric  surfaces  in  the  form 

Q(t)  =  HA^CTi4  -  T24) 

where 

Q  is  the  net  rate  of  radiation  heat  exchange 
H  is  the  radiation  factor 
A} ,A2  are  the  areas  of  the  surfaces 
T i ,T2  are  the  temperatures  of  the  surfaces. 

H  may  be  temperature  dependent  and  is  given  on  the  TABLEH1  card. 

The  CONVEC  card  defines  convective  heat  exchange  between  an  isopara¬ 
metric  surface  element  and  its  surroundings  in  the  form 
Q(t)  =  H,A(T  -  TJ 

where 

Q  is  the  heat  flow  per  unit  time 

H  is  the  film  coefficient 
c 

A  is  the  area  of  the  element 
T  is  the  temperature  of  the  surface 
T^  is  the  ambient  temperature. 
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Hc  may  be  time  dependent  and  is  given  on  the  TABLHC1  card.  may  be 
time  dependent  and  is  given  on  the  TABLET1  card.  (The  portion  of  A(t) 
defined  by  -H  ATn  is  the  load  added  to  the  load  vector.  The  portion 
defined  by  HcA  is  added  into  the  conduction  matrix.) 

The  HGEN  card  defines  internal  heat  generation  for  isoparametric 
elements  in  the  form 

Q(t)  =  GV 

where 

Q  is  the  load  to  be  applied 

G  is  the  heat  generation  per  unit  time  for  the  element 
V  is  the  volume  of  the  element. 

G  may  be  temperature  dependent  and  is  given  on  the  TABLEG1  card. 

The  HFLUX  card  defines  a  boundary  heat  input  per  unit  area  for  iso¬ 
parametric  elements  in  the  form 
Q(t)  -  qA 

where 

Q  is  the  load  to  be  applied 
q  is  the  boundary  heat  input  per  unit  area 
A  is  the  area  of  the  surface, 
q  may  be  time  dependent  and  is  given  on  the  TABLEQ1  card. 


1-5-7  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  APPISO  Isoparametric  Element  Problem  Type  Flag 

Description:  Defines  the  type  of  analysis  to  be  performed  using  iso¬ 
parametric  elements. 

Format  and  Example: 


L 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CODE 

SB 

0 

Field  Contents 


CODE 


Indicator 
of  type  of 
run 


0 

1 

2 


thermal-structural  combination 
structural  only 
thermal  statics  only 


Remarks  : 


1.  One  and  only  one  APPISO  bulk  data  card  must  be  present  when 
isoparametric  elements  are  used. 

2.  CODE=0  implies  that  Rigid  Format  14,  i.e.  thermal  transient- 
structural  statics  combination,  is  to  be  run.  Conduction 
and  capacitance  matrices  will  be  computed  in  the  thermal 
portion,  and  stiffness  matrices  will  be  computed  in  the 
structural  portion.  (If  only  a  thermal  transient  analysis 
is  desired,  CODE=0  must  be  specified.) 

3.  If  C0DE=1,  then  all  Rigid  Formats  except  4,  5,  6,  and  14 
(presently)  may  be  used  because  no  differential  stiffness 
matrices  nor  piecewise  linear  analysis  capabilities  have 
been  included  for  isoparametric  elements. 

4.  If  C0DE=2,  only  Rigid  Format  1  may  be  used. 


5.  An  APPISO  bulk  data  cord  must  be  present  in  thermal  problems 
with  BAR  elements. 
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BULK  DATA  DECK 


Input  Data  Card  CIS2D4  Quadrilateral  isoparametric  element  connection 

Description:  Defines  a  quadrilateral  isoparametric  membrane  element 
(IS2D4)  of  the  structural  model. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CIS2D4 

EID 

PID 

Gl 

G2 

G3 

G4 

ID1 

TH 

CIS2D4 

8 

2 

8 

6 

1 

10 

Field 


Contents 


EID  Element  identification  number  (Integer  >0) 

PID  Identification  number  of  a  PIS2D4 

property  card  (Integer  >0) 

Gl,  G2,  G3,  G4  Grid  point  identification  numbers  of 

connection  points  (Integers  >0;  Gl  thru 
G4  must  be  unique) 

ID1  Reserved  for  possible  later  use 

TH  Material  property  orientation  angle 

in  degrees  (Real)  The  sketch  below  gives 
the  sign  convention  for  TH 
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Remarks: 


.  Element  identification  numbers  must  be  unique  with  respect 
to  all  other  element  identification  numbers. 

.  Grid  points  G1  thru  G4  must  be  ordered  as  shown  above, 

.  This  element  is  a  planar  element,  i,e„  G1  thru  G4  must 
lie  in  a  plane. 

.  Stresses  are  computed  in  the  element  coordinate  system. 
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BULK  DATA  DECK 


Input  Data  Card  CIS2D8  Quadriparabolic  isoparametric  element  connection 

Description:  Defines  a  quadriparabolic  isoparametric  membrane  element 
(IS2D8)  of  the  structural  model. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CIS2D8 

BID 

PID 

Gl 

G2 

G3 

G4 

G5 

G6 

+abc 

CIS2D8 

16 

2 

12 

10 

15 

18 

22 

3 

+A 

+abc 

G7 

G8 

IDI 

TH 

+A 

7 

11 

Field 


Contents 


EID 


Element  identification  number  (Integer  >0) 


PID 


Identification  number  of  a  PIS2D8 
property  card  (Integer  >  0) 


Gl,  G2,  .  .  ,,G8  Grid  point  identification  numbers  of 

connection  points  (Integers  >0;  Gl 
thru  G8  must  be  unique) 

ID1  Number  of  Gauss  quadrature  points  (ID1=2  or 

3 — default  is  2). 
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TH 


Material  property  orientation  angle 
in  degrees  (Real) .  The  sketch  below 
gives  the  sign  convention  for  TH. 


Remarks : 


.  Element  identification  numbers  must  be  unique  with  respect 
to  all  other  element  identification  numbers. 

.  Grid  points  G1  thru  G8  must  be  ordered  as  shown  above. 

.  This  element  is  a  planar  element,  i.e.,  G1  thru  G8  must 
be  in  a  plane. 

.  Stresses  are  computed  in  the  element  coordinate  system. 


2 . 4-28 j  (  7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  CIS3D8  Solid  Isoparametric  Element  Connection 

Description:  Defines  a  solid,  8-grid-point  isoparametric  element 
(IS3D8)  of  the  structural  model  without  reference  to  a  property  card. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CIS3D8 

EID 

Gl 

G2 

G3 

G4 

G5 

G6 

G7 

+  abc 

CIS3D8 

5 

8 

9 

10 

11 

2 

3 

m 

+A 

+abc 

G8 

ID1 

ID2 

MID 

+A 

5 

1 

1 

Field  Contents 


EID  Element  identification  number  (Integer  >0) 

Gl,  G2,  G3,  .  .  ,G8  Grid  point  identification  numbers  of 

connection  points  (Integers  >0;  Gl 
thru  G8  must  be  unique) 

ID1  Coordinate  system  identification  number 

of  the  rectangular  coordinate  system 
which  defines  the  material  axes  (Integer 
>0;  0  means  the  basic  coordinate  system). 
(Not  required  nor  desired  for  isotropic 
materials.) 
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ID2 


Reserved  for  possible  later  use 


MID 


Remarks: 


Material  identification  number 
(Integer  > 0) 


.  Element  identification  numbers  must  be  unique  with  respect 
to  all  other  element  identification  numbers. 

.  Grid  points  G1  thru  G8  must  be  ordered  as  shown  above,  i.e. 
counterclockwise  on  opposite  face  as  viewed  from- face  1 
(G1-G4) .  Also,  G1  and  G5  must  share  the  same  edge. 

.  The  material  property  identification  number  must  reference 
only  a  MAT7  or  MAT9  card. 

.  ID1  must  refer  to  a  rectangular  coordinate  system. 

.  Stresses  are  computed  in  material  coordinate  system  ID1. 
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BULK  DATA  DECK 


Input  Bata  Card  CIS3D20  Solid  Isoparametric  Element  Connection 

Description:  Defines  a  solid,  20-grid-point  isoparametric  element 
(IS3D20)  of  the  structural  model  without  reference  to  a  property  card. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CIS3D20 

EID 

Gl 

G2 

G3 

G4 

G5 

G6 

G7 

+  abc 

CIS3D20 

6 

1 

2 

3 

4 

5 

6 

7 

+A 

+abc 

G8 

G9 

G10 

Gil 

G12 

G13 

G14 

G15 

+def 

+A 

8 

9 

10 

11 

12 

13 

14 

15 

+B 

+def 

G16 

G17 

G18 

G19 

G20 

ID1 

ID2 

MID 

+B 

16 

17 

18 

19 

20 

2 

3 

Field 


Contents 


EID  Element  identification  number  (Integer  >0) 

Gl,  G2,  .  .  ,,G20  Grid  point  identification  numbers 

of  connection  points  (Integers  >0;  Gl 
thru  G20  must  be  unique) 


2. 4-2 8m  (7/1/72-NSRDC) 


ID1 


ID2 

MID 


Remarks: 


Coordinate  system  identification 
number  of  the  rectangular  coordinate 
system  which  defines  the  material  axes 
(Integer  *0;  0  means  the  basic  coordinate 

system).  (Not  required  nor  desired  for 
isotropic  materials.) 

Number  of  Gauss  quadrature  points  (ID1=2  or  3- 
default  is  2) . 

Material  identification  number  (Integer  >0) 


.  Element  identification  numbers  must  be  unique  with  tespect 
to  all  other  element  identification  numbers. 

.  Grid  points  G1  thru  G20  must  be  ordered  as  shown  above. 

.  The  material  property  identification  number  must  reference 
only  a  MAT7  or  MAT9  card. 

.  ID1  must  refer  to  a  rectangular  coordinate  system. 

.  Stresses  are  computed  in  material  coordinate  system  ID1. 


2.4-28U  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  CONVEC  Convection  Specifications 

Description:  Specifies  convective  heat  exchange  between  a  surface  of 
an  isoparametric  element  and  its  surroundings. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CONVEC 

TH 

TT 

HC 

TC 

SI 

S2 

S3 

S4 

+  abc 

CONVEC 

3 

.05 

9 

10 

2 

6 

+A 

+abc 

S5 

S6 

S7 

S8 

S9 

S10 

Sll 

S12 

+def 

+A 

21 

18 

32 

ENDT 

(etc.) 

Field  Contents 


TH  Identification  number  of  a  TABLHC1 

card  (Integer  >0  or  blank) 

TT  Identification  number  of  a  TABLET1 

card  (Integer  >0  or  blank) 

HC  If  TH=0  or  blank,  HC  is  the  constant 

value  for  H  for  all  values  of  time  (Real) 


2.  4-42a  (7/1/72 -NSRDC) 


TC 


Si 


Remarks: 


If  TT=0  or  blank,  TC  is  the  constant 
value  for  the  ambient  temperature 
at  all  values  of  time  (Real) 

Identification  numbers  of  surface 
elements  or  two-dimensional  isopara¬ 
metric  elements  (Integers  >0) 

.  This  card  implies  the  relationship 

q  =  HcA(T-To) 

is  the  heat  flow  per  unit  time  due  to  convection 
exchange 

is  the  film  coefficient,  which  may  be  time  dependent 
is  the  area  of  the  surface 

is  the  temperature  of  the  surface,  which  is  given 
by  the  average  temperature  of  the  grid  points 
defining  the  surface. 

is  the  ambient  temperature,  which  may  be  time 
dependent 

.  The  term  H  A  can  be  considered  a  contribution  to  the  con- 
c 

duction  (stiffness)  matrix,  while  -HcATro  can  be  considered 

a  linear  external  load.  The  contribution  to  the  k.  .  term 

il 

of  the  conduction  matrix  for  element  Sk  is 
H  //N  N  dS 

c.  1  i 

and  the  linear  external  load  applied  to  grid  point  i  is 

H  T  J7N.  dS 
c  00  1 

where  N-TK  is  the  shape  function  of  the  surface  element 
applied  at  grid  points  i  and  j. 


where  Q 

H 

< 

A 

T 


Too 


2.4-42b  (7/1/72-NSRDC) 


3.  The  end  of  the  list  is  indicated  by  the  BCD  string 
"ENDT”  in  The  field  following  the  last  entry.  An  error 
is  detected  if  any  continuation  cards  follow  the  card 
containing  the  end-of-list  flag  "ENDT". 

4.  In  remark  2. ,  if  Hc  is  time-dependent,  the  program  will 
interrupt  the  integration  process  to  recompute  the 
contributions  to  the  conduction  matrix  at  each  time  step. 
This  could  be  a  time-consuming  process.  Therefore,  an 
option  is  included  whereby  no  contributions  are  made  to  the 
conduction  matrix,  but  the  load  applied  is  computed  to  be 

Q  =  H  A  (T  -  T  ) 

c  “ 

where  T  is  the  temperature  at  the  previous  time  step. 

While  a  much  less  time-consuming  process,  this  method 
could  cause  instabilities  in  the  integration  algorithm.  This 
option  is  implemented  by  including  a  PARAM  METCON  bulk 
data  card,  with  a  value  of  1,  in  the  bulk  data  deck. 


2.4-42c  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  HFLUX  Boundary  Heat  Input  Specification 

Description:  Specifies  boundary  heat  input  per  unit  area  on  a  surface 
of  an  isoparametric  element. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

HFLUX 

TQ 

QC 

SI 

S2 

S3 

S4 

S5 

S6 

+  abc 

HFLUX 

3 

9 

10 

2 

6 

21 

18 

+A 

+  abc 

S7 

S8 

S9 

S10 

Sll 

S12 

S13 

S14 

+def 

+A 

32 

ENDT 

(etc.) 


Field  Contents 

TQ  Identification  number  of  a  TABLEQ1 

card  (Integer  ^0  or  blank) 

QC  If  TQ=0  or  blank,  QC  is  the  constant 

value  of  q  for  all  time  values  (Real) 

Identification  numbers  of  surface 
elements  or  two-dimensional  isopara¬ 
metric  elements  (Integer  >0) 


2. 4-138g  (7/1/72 -NSRDC) 


Remarks :  1.  This  card  implies  the  relationship 


Q(t)  =  q  (t)A 

where  Q  is  the  total  load  applied  to  the  surface 
q  is  the  boundary  heat  input  per  unit  area 
A  is  the  area  of  the  surface 

2.  Q  will  be  treated  as  a  linear  load  applied  to  the  grid 

points  defining  the  surface.  The  load  applied  to  grid  point 
i  is 


q  (t)//N.ds 

where  N  is  the  shape  function  for  the  surface  applied 
i 

at  grid  point  i. 

3„  The  end  of  the  list  is  indicated  by  the  BCD  string  "ENDT" 
in  the  field  following  the  last  entry.  An  error  is 
detected  if  any  continuation  cards  follow  the  card  containing 
the  end-of-list  flag  "ENDT". 


2.  4-138h  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  HGEN  Internal  Heat  Generation  Specifications 

Description:  Specifies  internal  heat  generation  for  isoparametric 
elements . 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

HGEN 

TG 

GC 

EIDI 

EID2 

EID3 

EID4 

EID5 

EID6 

+abc 

HGEN 

3 

9 

10 

2 

6 

21 

18 

+A 

+abc 

EID7 

EID8 

EID9 

EID10 

EIDI  1 

EID12 

EID13 

EID14 

+def 

+A 

32 

ENDT 

(etc.) 


Field 


Contents 


TG 


Identification  number  of  a  TABLEG1 
card  (Integer  -0  or  blank) 


GC  If  TG=0  orblank  ,  GC  is  the  constant 

value  of  G  for  all  temperature  values 
(Real) 

EIDi  Element  identification  number  of  an 

isoparametric  element  (Integer >  0) 


2.4-1381  (7/1/72-NSRDC) 


Remarks:  1.  This  card  implies  the  relationship 


Q  =  GV 

where  Q  is  the  load  to  be  applied 

G  is  the  heat  generation  per  unit  time  for  the 
element  which  may  be  temperature  dependent 
V  is  the  volume  of  the  element 

2.  The  load  Q  to  be  applied  at  grid  point  i  of  an  element  is 


G(Tj)///N  dxdydz 

if  the  element  is  three-dimensional,  or 


G(Ti)h//Nidxdy 

if  the  element  is  two-dimensional, 

where  N  ^  is  the  shape  function  for  the  element  applied 
at  grid  point  i 

h  is  the  thickness  of  a  two-dimensional  element 
is  the  temperature  at  grid  point  i 

3.  The  end  of  the  list  is  indicated  by  the  BCD  string  "ENDT" 
in  the  field  following  the  last  entry.  An  error  is  detected 
if  any  continuation  cards  follow  the  card  containing  the 
end-of-list  flag  "ENDT". 


2.  4-138j  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MAT 7  Material  Property  Definition 

Description:  Defines  the  material  properties  for  linear,  temperature- 
independent,  isotropic  materials  for  isoparametric  elements. 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MAT  7 

MID 

E 

G 

NU 

RHO 

Q 

TREF 

GE 

+abc 

MAT  7 

16 

1 

0.3 

3.2 

mm 

125. 

.1 

+A 

+abc 

K 

C 

+A 

.31 

.60 

Field 


MID 

E 

G 

NU 

RHO 


Contents 

Material  identification  number  (Integer >  0) 
Young's  modules  (Real>0.0  or  blank) 

Shear  modulus  (Real >  0.0  or  blank) 

Poisson's  ratio  (Real  or  Blank) 

Mass  density  (Real) 


2.4-146e  (7/1/72-NSRDC) 


A 


Thermal  expansion  coefficient  (Real) 


TREF 

GE 

K 

C 

Remarks : 


Remarks 
(CODE  = 


Thermal  expansion  reference  temperature 
(Real) 

Structural  element  damping  coefficient 
(Real) 

Conductivity  coefficient  (Real) 

Specific  heat  (Real) 

1.  The  material  identification  number  must  be  unique  for  all 
MAT1,  MAT2,  MAT3,  MAT4,  MAT5,  MA.T7 ,  MAT8,  and  MAT  9  cards. 

2„  MAT 7  materials  may  be  made  temperature  dependent  by  use 
of  the  MATT 7  card. 

3.  The  continuation  card  is  not  required.  If  it  is  left  oui^ 

K  and  C  will  be  0.0. 

,  5,  and  6  apply  whenever  a  structural  analysis  is  to  be  run. 

or  1  on  the  APPISO  bulk  data  card.) 

4.  One  of  E  or  G  must  be  positive  (i.e.,  either  E  >0.0  or 
G  >0.0)  or  both  E  and  G  may  be  >0.0. 

5.  If  any  one  of  E,  G,  or  NU  is  blank,  it  will  be  computed  to 
satisfy  the  identity  E  =  2(1  +  NU)G;  otherwise,  values 
supplied  by  the  user  will  be  used. 

6.  The  mass  density  RHO  will  be  used  to  automatically  compute 
mass  in  a  structural  dynamics  problem.  In  a  thermal 
transient  problem,  RHO  *  C  will  be  used  as  the  coefficient 


2. 4-146f  (7/1/72-NSRDC) 


6.  (cont.) 

of  the  temperature  derivative  with  respect  to  time  in 
the  heat  flow  equation,  i.e.,  RH0*C  =  capacitance. 


2.  4-146g  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MATS  Material  Property  Definition 

Description:  Defines  the  material  properties  for  linear,  temperature- 
independent,  anisotropic  materials  for  two-dimensional  isoparametric 
elements . 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MAT8 

MID 

Gil 

G12 

G13 

G22 

G23 

G33 

RHO 

+abc 

MATS 

17 

6.2+3 

5.1+3 

3.2 

+A 

+abc 

Al 

A2 

A12 

TREF 

GE 

KX 

KY 

C 

+A 

.15 

125. 

.1 

.31 

.5 

.60 

Field  Contents 

MID  Material  identification  number  (Integer  >0) 

Gij  The  3x5  symmetric  material  property 

matrix  (Real) 

RHO  Mass  density  (Real) 


Ai 


Thermal  expansion  coefficient  vector 
(Real) 


TREF 


Thermal  expansion  reference  temperature 
(Real) 


2 . 4-146h  (7/ 1/72-NSRDC) 


GE 

KX,  KY 

C 

Remarks : 


Structural  element  damping  coefficient  (Real) 

Conductivity  coefficients  in  the 

x^and  y-directions,  respectively  (Real) 

Specific  heat  (Real) 

.  The  material  identification  numbers  must  be  unique  for  all 
MAT1,  MAT2,  MAT3,  MAT4,  MAT5,  MAT7,  MAT8,  and  MAT9  cards. 

.  MAT8  materials  may  be  made  temperature  dependent  by  use 
of  the  MATT8  card. 

.  The  mass  density  RHO  will  be  used  to  automatically  compute 
mass  in  a  structural  dynamics  problem.  In  a  thermal 
transient  problem,  RHO*C  will  be  used  as  the  coefficient 
of  the  temperature  derivative  with  respect  to  time  in 
the  heat  flow  equation,  i.e.,  RHO*C  •  capacitance. 


2.  4-146i  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MAT9  Material  Property  Definition 

Description;  Defines  the  material  properties  for  linear,  temperature- in¬ 
dependent  ,  anisotropic  materials  for  solid  isoparametric  elements. 

Format  and  Example: 


2.4-1463  (7/1/72-NSRDC) 


Field 


Contents 


MID 

Material  identification  number 
(Integer  >  0) 

Gij 

The  6X6  symmetric  material  prop¬ 
erty  matrix  (Real) 

RHO 

Mass  density  (Real) 

Ai 

Thermal  expansion  coefficient  vector 
(Real) 

TREF 

Thermal  expansion  reference  temperature 
(Real) 

KX,  KY,  KZ 

Conductivity  coefficients  in  the 

x- ,  y-,and  z-directions ,  respectively 

(Real) 

C 

Specific  heat  (Real) 

Remarks:  1.  The  material  identification  numbers  must  be  unique  for  all 
MAT1,  MAT 2,  MAT3,  MAT4,  MAT5,  MAT7 ,  MAT 8  and  MAT 9  cards. 

2.  MATg  materials  may  be  made  temperature  dependent  by  use  of 
the  MATT9  card, 

3.  The  mass  density  RHO  will  be  used  to  automatically  compute 
mass  in  a  structural  dynamics  problem.  In  a  thermal 
transient  problem,  RHO*C  will  be  used  as  the  coefficient 
of  the  temperature  derivative  with  respect  to  time  in  the 
heat  flow  equation,  i.e.,  RH0*C  =  capacitance. 


2.  4-l46k  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MATT 7  Material  Temperature  Dependence 

Description:  Specifies  table  references  for  material  properties  which 
are  temperature  dependent. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MATT  7 

MID 

Rl 

R2 

R3 

R4 

R5 

R6 

R7 

+abc 

MATT7 

16 

32 

18 

+A 

+abc 

R8 

R9 

+A 

12 

Field 


Contents 


MID  Material  property  identification 

number  which  matches  the  identifi¬ 
cation  number  on  some  basic  MAT7 
card  (Integer  >0) 

Ri  References  to  table  identification 

numbers  (Integers  -  0) 


2.4-  154a  (7/1/72-NSRDC) 


Remarks:  1.  Blank  or  zero  entries  mean  no  table  dependence  of  the 
referenced  quantity  on  the  basic  MAT7  card. 

2.  TABLEM1,  TABLEM2,  TABLEM3,  or  TABLEM4  type  tables  may  be 
used. 

3.  If  a  thermal  transient  analysis  (Rigid  Format  14)  is 
being  run,  the  appearance  on  a  MATT  7  card  of  a  non-zero 
Ri  for  the  K  and/or  C  variables  (R8,  R9)  must  be  accom¬ 
panied  by  a  TRANGE  bulk  data  card  for  that  Ri. 


2. 4-154 b  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MATT 8  Material  Temperature  Dependence 

Description:  Specifies  table  references  for  material  properties  which 
are  temperature  dependent. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MATT  8 

MID 

Rl 

R2 

R3 

R4 

R5 

R6 

R7 

+abc 

MATT  8 

17 

32 

15 

12 

+A 

+abc 

R8 

R9 

RIO 

Rl  1 

R12 

R13 

R14 

R15 

+A 

9 

2 

8 

Field 


Contents 


MID  Material  property  identification 

number  which  matches  the  identification 
number  on  some  basic  MAT 8  card 
(Integer  >0) 

Ri  References  to  table  identification 

number  (Integer  -0) 


2.4-154  c  (7/1/72-NSRDC) 


Remarks :  1.  Blank  or  zero  entries  mean  no  table  dependence  of  the 

referenced  quantity  on  the  basic  MAT8  card. 

2.  TABLEM1,  TABLEM2,  TABLEM3,  or  TABLEM4  type  tables  may  be 
used. 

3.  If  a  thermal  transient  analysis  (Rigid  Format  14)  is  being 
run,  the  appearance  on  a  MATT^  card  of  a  non- zero  Ri 

for  the  KX,  KY,  and/or  C  variables  (R13,  R14,  R15)  must 
be  accompanied  by  a  TRANGE  bulk  data  card  for  that  Ri. 


2.  4-154d  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  MATT 9  Material  Temperature  Dependence 

Description:  Specifies  table  references  for  material  properties  which 
are  temperature  dependent. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MATT  9 

MID 

R1 

R2 

R3 

R4 

R5 

R6 

R7 

+abc 

MATT9 

17 

32 

18 

17 

+A 

+abc 

R8 

R9 

RIO 

Rll 

R12 

R13 

R14 

R15 

+def 

+A 

12 

+  B 

+def 

R16 

R17 

R18 

R19 

R20 

R21 

R22 

R23 

+ghi 

+B 

5 

10 

+C 

+ghi 

R24 

R25 

R26 

R27 

R28 

R29 

R30 

R31 

+jkl 

+C 

9 

11 

+D 

+jkl 

R32 

R33 

R34 

+D 

8 

2 . 4- 154  e  (7/1/72-NSRDC) 


Field 


Contents 


MID 


Ri 


Remarks : 


Material  property  identification 
number  which  matches  the  identifi¬ 
cation  number  on  some  basic  MAT9 
card  (Integer  >0) 

References  to  table  identification 
number  (Integer  >0) 

.  Blank  or  zero  entries  mean  no  table  dependence  of  the 
referenced  quantity  on  the  basic  MAT9  card. 

.  TABLEM1,  TABLEM2,  TABLEM3,  TABLEM4  type  tables  may  be  used. 

.  If  a  thermal  transient  analysis  (Rigid  Format  14)  is  being 
run,  the  appearance  on  a  MATT9  aard  of  a  non- zero  Ri  for 
the  KX,  KY,  KZ,  and/or  C  variables  (R31  through  R34)  must 
be  accompanied  by  a  TRANGE  bulk  data  card  for  that  Ri. 


2.4-154f  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  PIS2D4  Quadrilateral  Isoparametric  Membrane  Property 

Description:  Used  to  define  the  properties  of  a  quadrilateral  isoparametric 
membrane.  Referenced  by  the  CIS2D4  card. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PIS2D4 

PID 

MID 

0 

PIS2D4 

2 

1 

0.5 

Field  Contents 


PID 


Property  identification  number 
(Integer  >  0) 


MID  Material  identification  number 

(Integer  >  0) 

T  Thickness  of  membrane  (Real) 


Remarks :  1.  All  PIS2D4  cards  must  have  unique  property  identification 

numbers . 

2.  The  material  property  identification  number  must  ref erence  only 
a  MAT7  or  MAT8  card. 


2.4-192a  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  PIS2D8  Quadriparabolic  Isoparametric  Membrane  Property 

Description:  Used  to  define  the  properties  of  a  quadriparabolic  isopara¬ 
metric  membrane.  Referenced  by  the  CIS2D8  card. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PIS2D8 

PID 

MID 

T 

|B 

PIS2D8 

2 

1 

0.5 

H 

Field  Contents 


PID 


Property  identification  number 
(Integer  >0) 


MID 


Material  identifisation  number  (Integer 
>0) 


T 


Thickness  of  membrane  (Real) 


Remarks :  1.  All  PIS2D8  cards  must  have  unique  property  identification 

members . 

2.  The  material  property  identification  number  must  reference  only 
a  MAT7  or  MAT8  card . 


2.  4- 192b  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  PL0AD2  Pressure  Load 

Description:  Defines  a  uniform  static  pressure  load  applied  to  two- 
dimensional  elements.  Only  IS2D4,  IS2D8,  QUAD1 ,  QUAD 2 , 

QDMEM ,  QDPLT ,  SHEAR,  SURF 4 ,  SURF8 ,  TRBSC ,  TRIA1,  TRIA2 ,  TRMEM,  TRPLT  or 
TWIST  elements  may  have  a  pressure  load  applied  to  them  via  this  card. 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PLOAD 2 

SID 

P 

EID 

EID 

EID 

EID 

EID 

EID 

PL0AD2 

21 

-3.6 

4 

16 

2 

Alternate  Form 


PL0AD2 

SID 

P 

EID1 

"THRU" 

EID2 

S3 

83 

SI 

PL0AD2 

i—l 

30.4 

16 

THRU 

48 

Field 


Contents 


SID 


Load  set  identification  number 
(Integer  >  0) 


P 


Pressure  value  (Real) 


EID  'i 
EID1  > 
EID2  ) 


Element  identification  number 
(Integer  >  0;  EID1  <  EID2) 


Remarks :  1.  EID  must  be  0  or  blank  for  omitted  entrys. 

2.  Load  sets  must  be  selected  in  the  Case  Control  Deck  (L0AD=SID) 
to  be  used  by  NASTRAN . 

3.  At  least  one  positive  EID  must  be  present  on  each  PL0AD2  card. 


4.  If  the  alternate  form  is  used,  all  elements  EID1  thru  EID2 
must  be  two-dimensional 


5.  The  pressure  load  is  computed  for  each  element  as  if  the 
grid  points  to  which  the  element  is  connected  were 
specified  on  a  PLOAD  card.  The  grid  point  sequence  specified 
on  the  element  connection  card  is  assumed  for  the  purpose 

of  computing  pressure  loads. 

6.  All  elements  referenced  must  exist. 


2.4-197  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  RADIAT1  Radiation  Specifications 

Description:  Specifies  radiation  heat  exchange  between  a  surface  of  an 
isoparametric  element  and  its  surroundings. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RAD I ATI 

TH 

TT 

HC 

TC 

SI 

S2 

S3 

S4 

+abc 

RAD TATI 

3 

50. 

9 

10 

2 

6 

+A 

+abc 

S5 

S6 

S7 

S8 

S9 

S10 

Sll 

S12 

+def 

+A 

21 

18 

32 

ENDT 

(etc.) 


Field  Contents 

TH  Identification  number  of  a 

TABLEH1  card  (Integer  -0  or  blank) 

TT  Identification  number  of  a  TABLET1 

card  (Integer  -0  or  blank  ) 

HC  If  TH=0  or  blank,  HC  is  the  constant 

value  of  H  for  all  temperature  values 
(Real) 


TC  If  TT  =  0  or  blank,  TC  is  the  constant 

value  for  the  ambient  temperature  at 
all  values  of  time  (Real) 


2.4-236 c  (7/1/72-NSRDC) 


Si 


Remarks : 


Identification  numbers  of  surface 
elements  or  two-dimensional 
isoparametric  elements.  (Integers  >0) 


.  This  card  implies  the  relationship 

Q  =  HACr4-^4) 


where  Q 

H 

A 

T 


T 

00 


is  the  net  rate  of  radiation  heat  exchange  between 
a  specified  surface  and  its  surroundings 
is  the  radiation  factor,  which  may  be  temperature 
dependent 

is  the  area  of  the  surface 

is  the  temperature  of  the  surface,  which  is  given 
by  the  average  temperature  of  the  grid  points 
defining  the  surface 

is  the  ambient  temperature,  which  may  be  time 
dependent 


.  Q  will  be  treated  as  a  non-linear  load  applied  to  the  grid 
points  defining  the  surface.  The  load  applied  to  grid  point  i 
is 


H(T.)CH-T  „4)//NidS, 


where  is  the  shape  function  for  the  surface  applied 
at  grid  point  i,  and  is  the  temperature  of  grid  point  i. 

The  end  of  the  list  is  indicated  by  the  BCD  string  "ENDT" 
in  the  field  following  the  last  entry.  An  error  is  detected 
if  any  continuation  cards  follow  the  card  containing  the 
end-of-list  flag  "ENDT". 


2.  4-236d  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  RAD I AT 2  Radiation  Specifications 

Description:  Specifies  radiation  heat  exchange  between  one  surface  of 
an  isoparametric  element  and  another  surface  of  an  isoparametric  element. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

RADIAT2 

TH 

HC 

SI 

S21 

S22 

S23 

S24 

S25 

+abc 

RAD I AT 2 

3 

9 

10 

2 

6 

21 

18 

+A 

+abc 

S26 

S27 

S28 

S29 

S210 

S211 

S212 

S213 

+def 

+A 

32 

ENDT 

(etc.) 

Field  Contents 


TH  Identification  number  of  a  TABLEH1 

card  (Integer  -0  or  blank) 

HC  If  TH  =  0  or  blank,  HC  is  the  constant 

value  of  the  radiation  factor  for 
all  temperatures  (Real) 


SI  Identification  number  of  a  surface 

element  or  two-dimBnsional  isoparametric 
element  from  which  the  radiation 
is  assumed  to  be  emanating  (Integer  >0) 


2.  4-236e  (7/1/72-NSRDC) 


S2i 


Remarks : 


Identification  number  of  a  surface 
element  or  two-dimensional 
isoparametric  element  to  which 
element  SI  is  assumed  to  be  radiating 
(Integer  >  0) 


.  This  card  implies  the  relationship 


Q  = 


»*lA2i 


4 

(Tj-T 


4 

2i) 


where  Q 

H 


is  the  net  rate  of  radiation  heat  exchange  be¬ 
tween  surface  SI  and  surface  S  2^ 
is  the  radiation  factor,  which  may  be  tem¬ 
perature  dependent 
is  the  area  of  SI 
is  the  area  of  S  2 

is  the  temperature  of  SI,  which  is  given  by 
the  average  temperature  of  the  grid  points 
defining  SI 

is  the  temperature  of  S2i,  which  is  given 
by  the  average  temperature  of  the  grid  points 
defining  S2^. 


Q  will  be  treated  as  a  non-linear  load  applied  to  the 
grid  points  defining  the  surfaces.  The  load  applied  to 
grid  point  j  of  surface  SI  is 

HCT^-T  4r)A2.//N.dS 

j  i  J 

where  Nj  is  the  shape  function  for  SI  applied  at  grid  point 
j.  The  load  applied  to  grid  point  j  of  surface  S^  is 

-HCT^-T^^A^NjdS 


2.4-236  f  (7/1/72-NSRDC) 


3.  The  roles  of  elements  SI  and  S2i  may  not  be  reversed 
on  any  other  RADIAT2  card. 

4.  The  end  of  the  list  is  indicated  by  the  BCD  string 
"ENDT"  in  the  field  following  the  last  entry.  An 
error  is  detected  if  any  continuation  cards  follow 
the  card  containing  the  end-of-list  flag  "ENDT". 


2.4-236g  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  SURF1  One-Dimensional  Surface  Element. 

Description:  Defines  a  one-dimensional  surface  element  of  the  structural 
model. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SURF1 

EID 

G1 

A 

SURF1 

3 

12 

oo 

Field  Contents 

EID  Element  identification  number 

(Integer  >  0) 


G1  Grid  point  identification  number 

(Integer  >  0) 

A  Area  associated  with  grid  point  G1 

(Real  >  0.0) 


Remarks :  1.  Element  identification  numbers  must  be  unique  with  respect 

to  all  other  element  identification  numbers. 

2.  Surface  elements  contribute  no  stiffness,  mass,  or  damping 
properties,  nor  will  any  stresses  or  forces  be  calculated 
for  them.  The  only  function  of  SURF1  surface  elements  is 
to  provide  a  convenient  way  of  specifying  radiation,  flux, 
and  convection  properties. 


2.4-270a  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  SURF4  Quadrilateral  Isoparametric  Surface  Connection 

Description:  Defines  a  quadrilateral  isoparametric  surface  element  of 
the  structural  model. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SURF4 

EID 

Gl 

G2 

G3 

G4 

SURF  4 

3 

6 

3 

1 

10 

Field  Contents 


EID 


Element  identification  number 
(Integer  >0) 


Gl,  G2,  G3,  G4  Grid  point  identification  numbers 

of  connection  points  (Integers  >0] 
Gl  thru  G4  must  be  unique) 


Remarks :  1.  Element  identification  numbers  must  be  unique  with  respect  to 

all  other  element  identification  numbers. 

2.  Grid  points  Gl  thru  G4  must  be  ordered  as  shown  in  the  sketch 
above . 


2.4-270*>  (7/1/72-NSRDC) 


3.  Surface  elements  contribute  no  stiffness,  mass,  or  damping 
properties,  nor  will  any  stresses  or  forces  be  calculated 
for  them.  The  only  function  of  SURF4  surface  elements  is 
to  provide  a  convenient  way  of  specifying  radiation,  flux,  and 
convection  properties  and  uniform  static  pressure  loads  on 
surfaces  of  isoparametric  elements. 


2.4-270c  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  SURF8  Quadriparabolic  Isoparametric  Surface  Element 

Description:  Defines  a  quadriparabolic  isoparametric  surface  element 
of  the  structural  model. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SURF  8 

EID 

Gl 

G2 

G3 

G4 

G5 

G6 

G7 

♦abc 

SURF  8 

3 

16 

22 

3 

19 

11 

8 

25 

+A 

+abc 

G8 

ID1 

+A 

18 

Field 


Contents 


EID  Element  identification  number 

(Integer  >0) 

Gl,  G2,  .  .  G8  Grid  point  identification  numbers 

of  connection  points  (Integers  >0; 

Gl  thru  G8  must  be  unique) 

ID1  Number  of  Gauss  quadrature  points  (ID1=2  or  3- 

default  is  2)  . 


2.4-270d  (7/1/72-NSRDC) 


Remarks:  1.  Element  identification  numbers  must  be  unique  with  respect 
to  all  other  element  identification  numbers. 

2.  Grid  points  G1  thru  G8  must  be  ordered  as  shown  in  the  above 
sketch. 

3.  Surface  elements  contribute  no  stiffness,  mass,  or  damping 
properties,  nor  will  any  stresses  or  forces  be  calculated 
for  them.  The  only  function  of  SURF8  surface  elements  is 

to  provide  a  convenient  way  of  specifying  radiation,  flux,  and 
convection  properties  and  uniform  static  pressure  loads  on 
surfaces  of  isoparametric  surfaces. 


2. 4-270  e  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TABLEG1  Internal  Heat  Generation  Factor  Specification 

Description:  Defines  a  tabular  function  for  temperature- dependent  internal 
heat  generation  factors.  Referenced  by  HGEN  bulk  data  cards. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

TABLEG1 

ID 

£21 

sa 

^1 

SSI 

ss 

+abc 

TABLEG1 

3 

+A 

+abc 

T1 

G1 

T2 

G2 

T3 

G3 

W 

G4 

+def 

+A 

0. 

.75 

100. 

.82 

200. 

.9 

ENDT 

(etc.) 


Field  Contents 

ID  Table  identification  number 

(Integer >  0) 

Ti,  Gi  Tabular  entries,  where  Ti  is  a 

temperature  value  (abscissa)  and 

Gi  is  the  internal  heat  generation 
factor  (ordinate)  (Real) 


2.4-280a  (7/1/72-NSRDC) 


Remarks:  1.  The  must  be  in  either  ascending  or  descending  order  but 
not  both. 

2.  Jumps  between  two  points  0\  =  T^+^)  are  allowed,  but  not 
at  the  end  points. 

3.  At  least  two  entries  must  be  present. 

4.  Any  T-G  entry  may  be  ignored  by  placing  the  BCD  string 
"SKIP"  in  either  of  the  two  fields  used  for  that  entry. 

5.  The  end  of  the  table  is  indicated  by  the  BCD  string  "ENDT" 
in  either  of  the  two  fields  following  the  last  entry. 

An  error  is  detected  if  any  continuation  cards  follow 
the  card  containing  the  end-of-table  flag  "ENDT". 

6.  At  temperature  T^,  the  value  of  G  is  G^.  Otherwise,  G=f(T), 
where  T  is  input  to  the  table  and  G  is  computed  by  linear 
interpolation  within  the  table  and  linear  extrapolation 
outside  the  table,  using  the  last  two  end  points  at  the 
appropriate  table  end.  At  jump  points,  the  average  is  used. 

7.  The  table  ID  must  be  unique  with  respect  to  all  TABLEH1, 
TABLHC1 ,  TABLET1 ,  TABLEQ1  and  TABLEG1  cards. 


2.4-280b  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TABLEH1  Radiation  Factor  Table 

Description:  Defines  a  tabular  function  for  temperature -dependent 
radiation  factors.  Referenced  by  RADIATI  and  RADIAT2  bulk  data  cards. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

TABLEH1 

ID 

X 

X 

X 

X 

X 

+abc 

TABLEH1 

3 

+A 

+abc 

T1 

Hi 

T2 

H2 

T3 

H3 

mm 

H4 

+def 

+A 

200 

.21 

850. 

.33 

1475. 

.61 

ENDT 

(etc.) 


Field  Contents 

ID  Table  identification  number 

(Integer  >0) 

T. ,  H.  Tabular  entries,  where  T.  is  a 

1’  1  1 

temperature  value  (abscissa)  and 
IT  is  the  corresponding  radiation 
factor  (ordinate)  (Real) 


2.4-280C  (7/1/72-NSRDC) 


Remarks :  1 . 


The  must  be  in  either  ascending  or  descending  order  but 
not  both. 

2.  Jumps  between  two  points  (J\  =  +^)  are  allowed,  but  not 

at  the  end  points. 

3.  At  least  two  entries  must  be  present. 

4.  Any  T-H  entry  may  be  ignored  by  placing  the  BCD  string 
"SKIP"  in  either  of  the  two  fields  used  for  that  entry. 

5.  The  end  of  the  table  is  indicated  by  the  BCD  string 
"ENDT"  in  either  of  the  two  fields  following  the  last 
entry.  An  error  is  detected  if  any  continuation  cards  follow 
the  card  containing  the  end-of-table  flag  "ENDT"- 

6.  At  temperature  T^,  the  value  of  H  is  H^.  Otherwise,  H=f(T), 
where  T  is  input  to  the  table  and  H  is  computed  by  linear 
interpolation  within  the  table  and  linear  extrapolation  out¬ 
side  the  table,  using  the  last  two  end  points  at  the  appro¬ 
priate  table  end.  At  jump  points,  the  average  is  used. 

7.  The  table  ID  must  be  unique  with  respect  to  all  TABLEH1,  TABLHC1, 
TABLET1 ,  TABLEQ1  and  TAB LEG 1  cards. 


2. 4-280d  (7/1/72 -NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TABLEQ1  Boundary  Heat  Input  Value  Specification 

Description:  Defines  a  tabular  function  for  the  time-dependent  boundary 
heat  input.  Referenced  by  HFLUX  bulk  data  cards. 


(etc.) 


Field 


Contents 


Table  identification  number  (Integer  >  0) 

t^,  Tabular  entries,  where  t^  is  a  time 

value  (abscissa)  and  is  the  boundary 
heat  input  per  unit  area  (ordinate)  (Real) 


2.  4-288a  (7/1/72 -NSRDC) 


Remarks:  1.  The  t  ^  must  be  in  either  ascending  or  descending 

order  but  not  both. 

2.  Jumps  between  two  points  (t^  =  t  ^ )  are  allowed, 
but  not  at  the  end  points. 

3.  At  least  two  entries  must  be  present. 

4.  Any  t-q  entry  may  be  ignored  by  placing  the  BCD  string 
"SKIP"  in  either  of  the  two  fields  used  for  that 
entry. 

5.  The  end  of  the  table  is  indicated  by  the  BCD  string 
"ENDT"  in  either  of  the  two  fields  following  the 
last  entry.  An  error  is  detected  if  any  continuation 
cards  follow  the  card  containing  the  end-of-table 
flag  "ENDT". 

6.  At  temperature  t^,  the  value  of  q  is  q^.  Otherwise, 
q  =  f(t),  where  t  is  input  to  the  table  and  q  is 
computed  by  linear  interpolation  within  the  table  and 
linear  extrapolation  outside  the  table  using  the  last 
two  end  points  at  the  appropriate  table  end.  At 
jump  points,  the  average  is  used. 

7.  The  table  ID  must  be  unique  with  respect  to  all 
TABLEH1 ,  TABLHC1,  TABLET1 ,  TABLEQ1 ,  and  TABLEG1  cards. 


2.4-288b  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TABLET1  Ambient  Temperature  Table 

Description;  Defines  a  tabular  function  for  time-dependent  ambient 
temperatures  used  in  convection  and  radiation  calculations.  Re¬ 
ferenced  by  CONVEC  and  RADIAT1  bulk  data  cards. 

Format  and  Example: 


123456789  10 


T 

CO  1 

*2 

CM 

H8 

t3 

T 

00  3 

B 

Q 

+def 

■ 

0. 

1200. 

.5 

1325. 

1.25 

1950. 

ENDT 

(etc.) 


Field 

ID 


Contents 

Table  identification  number  (Integer  >0) 

Tabular  entries,  where  t  is  a  time 

i 

value  (abscissa)  and  To,  ^  is  the  value  of 
the  ambient  temperature  (ordinate)  (Real) 


2.4-290a  (7/1/72-NSRDC) 


Remarks: 


1.  The  must  be  in  either  ascending  or  descending  order 
but  not  both, 

2.  Jumps  between  two  points  (t^=  t^  +  are  allowed, 
but  not  at  the  end  points, 

3.  At  least  two  entries  must  be  present. 

4.  Any  t-T^  entry  may  be  ignored  by  placing  the  BCD 
string  "SKIP"  in  either  of  the  two  fields  used  for 
that  entry. 

5.  The  end  of  the  table  is  indicated  by  the  BCD  string 
"ENDT"  in  either  of  the  two  fields  following  the  last 
entry.  An  error  is  detected  if  any  continuation  cards 
follow  the  card  containing  the  end-of-table  flag  "ENDT" 

6.  At  time  t. ,  the  value  of  T  is  T*,  . .  Otherwise,  T^  = 

l  i 

f(t),  where  t  is  input  to  the  table  and  T^  is  computed 
by  linear  interpolation  within  the  table  and  linear 
extrapolation  outside  the  table,  using  the  last  two  end 
points  at  the  appropriate  table  end.  At  jump  points, 
the  average  is  used. 

7.  The  table  ID  must  be  unique  with  respect  to  all  TABLEH1 
TABLHC1.  TABLET1 ,  TABLEQ1 ,  and  TABLEG1  cards. 


2.4-290b  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TABLHC1  Convection  Film  Coefficient  Speci¬ 
fication 

Description:  Defines  a  tabular  function  for  time-dependent 
convection  film  coefficients.  Referenced  by  CONVEC  bulk  data 
cards . 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

IIIIII2I 

ID 

S3 

S3 

S3! 

S3 

S3 

^1 

83 

+abc 

SMlIl 

jiisiiiissi 

3 

+A 

+abc 

tl 

H 

C  1 

t2 

H 

c2 

t 

3 

H 

C3 

B 

H 

C4 

+def 

+A 

0. 

.2 

25. 

.33 

60. 

.5 

ENDT 

(etc.) 


Contents 

Table  identification  number  (Inte¬ 
ger  >  0) 

Tabular  entries,  where  t^  is  a 
time  value  >0; abscissa  and  H 

c . 

is  the  corresponding  convection 
film  coefficient  (ordinate)  (Real) 


Field 

ID 


2.4-290c  (7/1/72-NSRDC) 


Remarks :  1 .  The  must  be  in  either  ascending  or  descending 
order  but  not  both. 


2.  Jumps  between  two  points  (t^=  t^)  are  allowed, 
but  not  at  the  end  points. 

At  least  two  entries  must  be  present. 

4.  Any  t-Hc  entry  may  be  ignored  by  placing  the  BCD  string 
"SKIP"  in  either  of  the  two  fields  used  for  that  entry. 

5.  The  end  of  the  table  is  indicated  by  the  BCD  string  "ENDT" 
in  either  of  the  two  fields  following  the  last  entry.  An 
error  is  detected  if  any  continuation  cards  follow  the 
card  containing  the  end-of-table  flag  "ENDT". 

6.  At  temperature  t  the  value  of  Hc  is  Hc  .  Otherwise, 
Hc=f(t),  where  t  is  input  to  the  table  1  and  Hc  is 
computed  by  linear  interpo  lation  within  the  table 
and  linear  extrapolation  outside  the  table,  using  the 
last  two  end  points  at  the  appropriate  table  end.  At 
jump  points,  the  average  is  used. 

7.  The  table  ID  must  be  unique  with  respect  to  all 
TABLEH1,  TABLHC1 ,  TABLET1 ,  TABLEQ1,  and  TABLEG1  cards. 


2 . 4-290d  (7/1/72-NSRDC) 


BULK  DATA  DECK 


Input  Data  Card  TRANGE  Temperature  Range  Specifications 

Description:  Specifies  temperature  ranges  within  which  temperature- 
dependent  thermal  material  properties  are  to  remain  constant  during  a 
thermal  transient  analysis. 

Format  and  Example: 


123456789  10 


2m 

ID 

AT 

B 

T 

2 

T3 

B 

T5 

T, 

6 

+abc 

■g 

8 

5. 

100. 

200. 

300. 

400. 

500. 

600. 

+A 

(etc.) 


Contents 

Table  identification  number  (In¬ 
teger  >0) 

Tolerance  (Real  *  0.0) 

Temperature  values  (Real) 


Field 

ID 

AT 

T. 


2.4-306a  (7/1/72-NSRDC) 


Remarks:  1.  TRANGE  cards  will  be  used  during  a  thermal  transient 
analysis  to  determine  whether  a  change  in  element 
material  properties  is  required. 

2.  The  table  identification  number  ID  must  match  the 

table  identification  number  on  a  TABLEM1,  TABLEM2,  TABLEM3, 
or  TABLEM4  bulk  data  card  which  defines  a  tabular  function 
for  a  thermal  material  property  KX,  KY,  KZ,  or  C  on  a 
MAT7,  MAT8,  or  MAT9  bulk  data  card. 


3.  If  a  thermal  material  property  KX,  KY,  KZ,  or  C  is  defined 
to  be  temperature  dependent  (by  a  non- zero  entry  in 

the  appropriate  field  on  a  MATT7,  MATT8,  or  MATT9  bulk 
data  card) ,  then  a  TRANGE  card  specifying  the  appropriate 
table  identification  number  must  be  present. 

4.  The  T^  must  be  specified  in  strictly  increasing  order. 

5.  At  least  two  T^'s  must  be  specified. 

6.  The  end  of  the  list  is  indicated  by  the  BCD  string  "ENDT" 
in  the  field  following  the  last  entry.  An  error  is 
detected  if  any  continuation  cards  follow  the  card  con¬ 
taining  the  end-of-list  flag  "ENDT". 

7.  If  the  temperature  of  an  element  (as  given  by  the  average 
temperature  of  its  grid  points)  is  in  the  interval 

(Tj  +  AT,  T^  -  !  -  AT),  then  the  temperature  range  will 
be  considered  to  be  [T^,  T^  +  ^] ,  If  the  temperature 
of  the  element  at  the  previous  time  step  in  the  transient 
analysis  was  not  in  this  range,  new  material  properties 
for  the  element  will  be  computed.  Otherwise,  the  same 
properties  will  be  used. 
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8. 


If  the  temperature  of  an  element  is  in  the  range 
[T\  -  AT,  +  AT],  then  a  determination  will  be 
made  by  NASTRAN  as  to  whether  new  material  pro¬ 
perties  for  the  element  will  be  computed.  The 
program  will  attempt  to  minimize  the  number  of 
material  property  changes.  Therefore,  the  tolerance 
AT  gives  NASTRAN  further  leeway  in  determining  material 
property  changes.  AT  may  be  0.  However,  the  tempera¬ 
ture  ranges  may  not  be  null  intervals. 

9.  When  NASTRAN  determines  that  the  temperature  of  an 
element  lies  in  a  particular  interval,  the  value  of 
the  appropriate  material  property  is  the  average  of 
the  values  of  the  material  property  at  the  end  points 
of  the  range. 

10.  Temperatures  less  than  T^  are  assumed  to  lie  in  in¬ 
terval  [T^,  T2]  while  temperatures  greater  than  TN 
are  assumed  to  lie  in  the  interval  [T^  ^ ,  TN] 
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RIGID  FORMATS 


restarts.  The  deletion  of  operations  for  each  subset  is  controlled  by  the 
restart  tables. 

If  the  user  wishes  to  modify  the  DMAP  sequence  of  a  rigid  format 
in  some  manner  not  provided  for  in  the  available  subsets,  he  can  use  the 
ALTER  feature  described  in  Section  2.  Typical  uses  are  to  schedule  an 
EXIT  prior  to  completion,  in  order  to  check  intermediate  output, 
schedule  the  printing  of  a  table  or  matrix  for  a  diganostic  purposes, 
and  to  delete,  or  add  a  functional  module  to  the  DMAP  sequence.  Any 
DMAP  instructions  that  are  added  to  a  rigid  format  are  automatically 
executed  when  a  restart  is  performed.  The  user  should  be  familiar 
with  the  rules  for  DMAP  programming,  as  described  in  Section  5,  prior  to 
making  alterations  to  a  rigid  format. 

The  following  rigid  formats  are  currently  included  in  NASTRAN: 

1.  Status  Analysis 

2.  Static  Analysis  with  Inertia  Relief 

3.  Normal  Mode  Analysis 

4.  Static  Analysis  with  Differential  Stiffness 

5 .  Buckling  Analysis 

6.  Piecewise  Linear  Analysis 

7.  Direct  Complex  Eigenvalue  Analysis 

8.  Direct  Frequency  and  Random  Response 

9.  Direct  Transient  Response 

10.  Modal  Complex  Eigenvalue  Analysis 

11.  Modal  Frequency  and  Random  Response 

12 .  Modal  Transient  Response 

13.  Normal  Modes  with  Differential  Stiffness 

14.  Thermal  Transient  -  Structural  Static  Combination 

3.1.1  Input  File  Process 

The  Input  File  Processor  operates  in  the  Preface  prior  to  the 
execution  of  the  DMAP  operations  in  the  rigid  format.  A  complete 
description  of  the  operations  in  the  Preface  is  given  in  the  Programmer's 
Manual.  The  main  interest  here  is  to  indicate  the  source  of  data  blocks 
that  are  created  in  the  Preface  and  hence  appear  only  as  inputs  in  the 
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DMAP  sequences  of  the  rigid  formats.  None  of  the  data  blocks  created 
by  the  Input  File  Processor  are  checkpointed,  as  they  are  always 
regenerated  on  restart.  The  Input  File  Processor  is  divided  into  four 
parts.  The  first  part  (IFP1)  processes  the  Case  Control  Deck,  the  second 
part  (IFP)  processes  the  Bulk  Data  Deck,  the  third  part  (IFP3)  performs 
additional  processing  of  the  bulk  data  cards  associated  with  the  conical 
shell  element,  the  fourth  part  (IFP4)  performs  additional  processing 
of  the  bulk  data. 
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GENERAL  DESCRIPTION  OF  RIGID  FORMATS 


2.  WTMASS  -  optional  in  all  rigid  formats.  The  terms  of  the  structural 
mass  matrix  are  multiplied  by  the  real  value  of  this  parameter 
when  they  are  generated  in  SMA2.  Not  recommended  for  use  in 
hydroelastic  problems. 

3.  IRES  -  optional  in  all  statics  problems  (rigid  formats  1,  2,  4,  5 
and  6).  A  positive  integer  value  of  this  parameter  will  cause  the 
printing  of  the  residual  vectors  following  each  execution  of  SSG3. 

4.  LFREQ  and  HFREQ  -  required  in  all  modal  formulations  of  dynamics 
problems  (rigid  formats  10,  11  and  12)  unless  LMODES  is  used.  The 
real  values  of  these  parameters  give  the  frequency  range  (LFREQ  is 
lower  limit  and  HFREQ  is  upper  limit)  of  the  modes  to  be  used  in 
the  modal  formulation. 

5.  LMODES  -  required  in  all  modal  formulations  of  dynamic  problems 
(rigid  formats  10,  11  and  12)  unless  LFREQ  and  HFREQ  are  used. 

The  integer  value  of  this  parameter  is  the  number  of  lowest  modes 
to  be  used  in  the  modal  formulation. 

6.  G  -  optional  in  the  direct  formulation  of  all  dynamics  problems 
(rigid  formats  7,  8  and  9).  The  real  value  of  this  parameter  is 
used  as  a  uniform  structural  damping  coefficient  in  the  direct 
formulation  of  dynamics  problems.  Not  recommended  for  use  in 
hydroelastic  problems. 

7.  W3  and  W4  -  optional  in  the  direct  formulation  of  transient  response 
problems.  The  real  values  of  these  parameters  are  used  as  pivotal 
frequencies  for  uniform  structural  damping  and  element  structural 
damping  respectively.  The  parameter  W3  should  not  be  used  for 
hydroelastic  problems. 

8.  MOD ACC  -  optional  in  the  modal  formulation  of  frequency  response  and 
transient  response  problems.  A  positive  integer  value  of  this 
parameter  causes  the  Dynamic  Data  Recovery  module  to  use  the  mode 
acceleration  method.  Not  recommended  for  use  in  hydroelastic 
problems . 

9.  COUPMASS  -  optional  in  all  rigid  formats.  A  positive  integer  value 
of  this  parameter  will  cause  the  generation  of  coupled  mass  matrices 
rather  than  lumped  mass  matrices  for  all  bar  elements,  rod  elements, 
and  plate  elements  that  include  bending  stiffness.  This  option 
applies  to  both  structural  and  nonstructural  mass  for  the  following 
elements:  BAR,  CONROD,  QUAD1,  QUAD 2 ,  ROD,  TRIA1 ,  TRIA2 ,  TUBE.  Also, 
for  IS2D4,  IS2D8,  IS3D8,  and  IS3D20  elements  coupled  mass  and 
capacitance  matrices  will  be  computed.  In  Rigid  Format  14  coupled 
capacitance  matrices  will  be  computed  for  BAR  elements.  Since 
structural  mass  is  not  defined  for  the  following  list  of  elements, 

the  option  applies  only  to  the  nonstructural  mass:  QDPLT,  TRBSC,  TRPLT. 
A  negative  value  causes  the  generation  of  lumped  mass  and  capacitance 
matrices  for  all  the  above  elements.  (This  is  the  default.)  A  zero 
value  activates  the  following  parameters  described  under  10. 
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10.  CP BAR,  CPROD,  CPQUADl ,  CPQUAD2 ,  CPTRIA1,  CPTRIA2 ,  CPTUBE,  CPQDPLT , 
CPTRPLT,  CPTRBSC — optional  in  all  rigid  formats.  These  parameters 
are  active  only  is  COUPMASS=0.  A  positive  value  will  cause  the 
generation  of  coupled  mass  matrices  for  all  elements  of  that 
particular  type  as  shown  by  the  following  table: 


Parameter 


Element  Types 


CP  BAR 

CPROD 

CPQUADl 

CPQUAD2 

CPTRIA1 

CPTRIA2 

CPTUBE 

CPQDPLT 

CPTRBSC 


BAR 

ROD,  C0NR0D 

QUAD1 

QUAD2 

TRIA1 

TRIA2 

TUBR 

QDPLT 

TRBSC 


A  negative  value  (the  default)  for  these  parameters  will  cause  the 
general  of  the  lumped  mass  matrices  for  these  element  types. 


11.  DEC0M0PT  -  optional  for  frequency  response  problems.  The  integer 
value  of  this  parameter  is  used  to  control  the  type  of  arithmetic 
used  in  the  decomposition  of  the  dynamic  equations.  A  value  of  1 
(default)  means  that  double  precision,  complex  arithmetic  with  partial 
pivoting  will  be  used.  A  value  of  2  means  that  double  precision, 
complex  arithmetic  without  pivoting  will  be  used.  A  value  of  4 
means  that  single  precision,  complex  arithmetic  without  pivoting 
will  be  used. 


12.  OPTION  -  optional  for  static  analysis.  The  value  HEAT  is  used  to 
select  the  heat  transfer  option  for  Rigid  Format  No.  1. 

13.  BT  -  required  in  thermal  transient  problems  in  which  radiation  is 
being  used  and  temperatures  are  not  specified  in  °R  or  °K.  The  usual 
values  for  BT  are  459.69  and  273.16  for  specification  in  °F  and  °C, 
respectively. 

14.  METCQN  -  Optional  for  Rigid  Format  14.  A  value  of  1  will  cause 
the  program  to  make  convection  contributions  only  to  the  load  vector. 
For  full  details  see  Remark  4,  page  2.  4- 42c. 
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3. 15  THERMAL  TRANSIENT-STRUCTURAL  STATICS  COMBINATION 

3. 15. 1  DMAP  Sequence  for  Thermal  Transient-Structural  Statics  Combination 

RIGID  FORMAT  14 

NAS  TRAM  SOURCE  PROGRAM  COMPILATION 
DMAP-DMAP  INSTRUCTION 
NO. 

1  BEGIN  NO. 14  THERMAL  TRANSIENT -STRUCTURAL  STATICS  -  SERIES  M  S 

2  FILE  LLL«TAPE  t 

3  FILE  Q6*  A PPENO/PGG* APPENDAJGV* APPEND/ 6M“SAVE/KNN«SAVE  S 

4  PARAM  //C,N,NOP/V,N,INl«-l  S 

5  PARAM  //C.N.NOP/V.N.IPl*!  » 

6  PARAM  //C, N,NOP/V,N,ICASF*2  S 

7  MODA  /KIO  .GHO.KGGxy.GOO/C.N.O.O/C.N.O.O/C.N.O.O/C.N.O.O/C.N, 0.0/C.N, 

0/C. N,0/C,N,0/C,N,0/C,N,0/C,N,0 . 0/V, N, IFIR*0/V ,N ,IFIRST*0  S 

8  SAVE  IFIR.IFIRST  $ 

9  MOOA  /CASEXZ.t ./C.N.O . OVC, N /0.0/C.N, 0.0/C, N, 0.0/C, N.Q.O/C.N, 0/C, N.O/ 

C.N, O/C.N.O/C.N.O /D.N, 0 .O/V.N.IFIR/V.N.IFIRST  S 

10  PARAM  //C. N.NOP/V »  N, IF0URN4S 

11  PURGE  MGS, MNN.MFP.MAA/IMl  » 

12  PURGE  SKSGX.SGPST , SMG6 « SOGPWQ, SKGG .SRG.SYS, SUSET , SOGPST.SGM.SKNN, 

SKPf  » SKFS .SKSS.SGO.SKAA.SKOOB.SLOO.SUOO.SKLL.SKLR.SKRR.SLLL, 
SULL.SDM,  SPG,  SQR,  SPOT'S  RS.SPL. SUL  V.SUOOV.SRULV,  SR  UOV.SUGV.SPGG, 
SOS, OPGl,OQGl,OUGVt,OE9Sl,OEFFl,PUGVl,OESAVG,PLTX2, CASE XY, 

GEOH  3X» ESTZ  «  ECPTZ.SPTTX/ INI  $ 

13  GP1  GE8M1,GE0^2,/GPL,EQ€XIN,GPDT,CSTM,BGPOT,SIL/V,N,LUSET/  C,N%12V 

V,  N , NOGPOT  S 

14  SAVE  LUS€T , NOS PD TO 

15  PURGE  USET.GM,G0,KAA,BAA^MA1AVK4AA.PST,KFS.QP,EST/N0GPDT  S 

16  CHKPNT  GPL  <  EQEXIN.GPDT .CSTM.BGPOT.SIL, USET , GM.GO.K AA, BA A. MA A .K4 AA, 

PST , KFS.QP.EST  S 

17  COND  LBL 5, NOGPOT t 

18  GP2  GEBM2.EQFXIN/ECT  $ 

19  CHKPNT  ECT  S 

20  PLT9ET  PDOB  .EQEXIN.ECT/PLTtSETX  ,  PLTPAR.  GPSETS,  ELSETS/V.N.NSIL/  V.N, 

JUNPPLOT  «  -1  * 
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RIGID  FORMAT  14 

N  A  S  T  R  A  M  SOURCE  PM!  O  G  R  AM  COMPILATION 
OMAP-OMAP  INSTRUCTION 

NO. 

21 

SAVE 

NSIL . JUMPPLOT  $ 

22 

PRTMSG 

PLT6ETX//S 

23 

CHKPNT 

PLTPAR,GPSE7S*ELSETS  t 

24 

SETVAL 

//V,  N,PLTFLG/C,N,1/1V,N,PFILE/C,N,0  S 

25 

SAVE 

PLTFLGfPFILE  S 

26 

COND 

PI,  JUMP  PLOTS 

27 

PLOT 

PLT  PAR,GPSETS,ELSETfS,C*SECC,  BGPD T ,EGEX IN,SI L » , /PLOT  XI /  V,N, 
NSIL/V,N,iLUSET/V,N,UUMPPL0T/V,N,PLTFLG/V,N,  PFILE  s 

28 

SAVE 

JUMPPLOT, PLTFLG, PFILE  S 

29 

PRTMSG 

PLOT  X1//S 

30 

LABEL 

PI  S 

31 

GP3 

GEOM  3,EQE'X  IN  ,  GE0N2,'BGPDT/SLT,6PTT/C,N»123/V»N,N0GRAV/C»N,123  S 

32 

SAVE 

NOS  R  AV  S 

13 

PARAM 

//C',N,ANO/V,  N«  SKPNGG/V  VN*NOGRAV/V»Y»GRDPNT  $ 

34 

PURGE 

SMSG/SKPMGG  S 

35 

CHKPNT 

GPTfT  ,MGG,SMGG,SLT  $ 

36 

TA1, 

,  ECT  ,EPT  ,BGPDT  ,SIL»tGPT?»CSTM/EST,,GEI*ECPT,GPCT/V*N,LUSET/  C  ,  N, 
123/V,N,N09IMP*-1/C,N,D/V,N, NOGE NL*-1/V, N, GENEL  S 

37 

SAVE 

NOStI  HP,  NOGENL  ,  GENEUS 

38 

PARAM 

//C  ,  N, ANO/V ,  N, NOELHT /WN»NOGENL/V,N,NOSIMP  $ 

39 

CONO 

ERR‘0R6,  NOELHT  S 

40 

PURGE 

K4SG ,GPST ,OGPST , MGG  *BGG,  K4NN, K4FF, K4AA, MNN, MFF, MA A ,BNN , BFF , 
BAA, KGGX/NOSIHP/  OG PST /G ENEL S 

41 

CHKPNT 

EST , ECPT, GPCT»GEI*K4GG<6 PST  *  MGG, BGG, KGGX,OGPST ,  K4NN, K4FF ,K4 A  A, 
MNN  ,  MFF  ,MAA  ,  8NN,  BFFI,B*AA  S 

42 

DPDA 

Olf, MPT/O ITC.MATTArV^M, MAT  S 

43 

SAVE 

MAT'  S 
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RIGID  FORMAT  14 

NASTRAW  SOURCE  PI  H  G  R  A  H  COMPILATION 
OMAP-OMAP  INSTRUCTION 


NO. 

44 

PARAM 

//C,  N.NOT/V »  N,  NOMAT./V,  N .MAT  $ 

45 

EQUIV 

OITC  .OIT/NOMAT  S 

46 

CHKF'NT 

MATT  A  $ 

47 

CONO 

LBL110.IM1  $ 

48 

TABPT 

OITC,.,.//  S 

49 

LABEL 

LBL'110  t 

50 

INITEM 

CAS€CC,OYNAMICS,ESTi,ECRT  ,SIL  ,MATTA,DIT » EQEX IN/ES TX , ECPTX , 

MATT  AB/V, N,S CAL AR/V.N, MAT/ V,N, STEP/ V.N.NOTIC  $ 

51 

SAVE 

SCALAR, MAT, STEP, NOTUC  S 

52 

EQUIV 

ESfiX.EST/NOT  IC /ECPTX  .ECPT/NOTIC  t 

53 

CHKPNT 

ESTi,  ECPTVNATTAB  S 

54 

CONO 

LBL811.IH1  S 

55 

TABPT 

ESfOC  ,ECPTX,  ,  ,//  t 

56 

LABEL 

LBL811  t 

57 

PURGE 

MATT  A/IM1  $ 

58 

CONO 

LABL 82, IFIR  $ 

59 

JUMP 

LABL 83  $ 

60 

LABEL 

LABL  83  % 

64 

COND 

LABL 85, KC ON  S 

62 

LABEL 

LABL  82  % 

63 

CONO 

LBL1 ,NOSINP$ 

64 

SMAL 

CSTM ,NPT,eCPT,CPCT^  DIP/KGGX , K4GG ,GPST/V , N, NOGENL/  V, N, N0K4GGS 

65 

SAVE- 

N0K4GGS 

66 

PUROE 

K4NN  ,K4GG,K4EF,K4AA/N0K4GG$ 

6.7 

PUROE 

K4661,K4GG2 , K4GG3.K496N/N0K4GG  $ 

6« 

CHKPNT 

KGf  X ,GPST , K4GG, K4NW.K4PF , K4 A  A  $ 
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RIGID  FORMAT  14 

N  A  S  T  R  A  *  SOURCE  Pi  R  0  G  R  A  M  COMPILATION 

DHAP-OMAP  INSTRUCTION 

NO. 

59 

LABEL 

LABL85  % 

TO 

CONO 

LABL  90  *  IFIIR  $ 

71 

CONI? 

LABL  98,  KCAP  $ 

72 

LABEL 

LABL90  S 

73 

SMA2 

cstm,mpt,ecpt,gpct,oit/ngg,bgg/v,y,mtmass«i.o/v, N.N0MGG7  v,n, 
NOBGG*-i/V,  Y,COUPMAiS9=-i  S 

74 

SAVE 

NOItGG.NOBGGJ 

75 

PURGE 

BN»,  BFF  ,8AA ),  BGG/NOB,GQ* 

76 

JUMP 

LABL  98  S 

77 

LABEL 

LABL  98  $ 

78 

EQUI  V 

KGSX  ,KGGXY/IM1  t 

79 

CHKPNT 

MGS  ,  BGG  , BNN  ,  BFF  ,  BAA.',  MNN  ,  MFF  ,  MAA  S 

SO 

LABEL 

LBL1  S 

81 

EQUI  V 

KGSX  ,  KGG7NOGENL  $ 

82 

CHKPNT 

KGS  S 

83 

CONO 

LBlll.NOGENL  $ 

84 

SMA3 

GEII,  KGGX/KGG/ V,  N,  LUSET  7  V  ,  N,  NOGENL/V  ,N,  NOS  IMP  $ 

85 

CHKPNT 

KGS  S 

86 

LABEL 

LBL11  I 

87 

EQUI  V 

KGSX Y«KG£X7 1 PI  S 

88 

PARAH 

77~, N,MPY/y,N,NSKIPyO,N,0/C,N,0  S 

89 

COND 

LABL  91, IFIRST  S 

90 

GP4 

CASECC, GE0M4 ,EQE  XIM, SIL , GPDT/RG, YS,USET,/V,N,LUS ET/V.N, MPCFi* 
-17  V  ,  N,  MPCF2“-17 VjNUSINGLE*- 1/V  ,N*OMIT*-l/V,N»REACT»— 1/V»N, 
NSKIP/V,N,REPEAT/V,INWN0SET=-17V,N,N0L/V,N,N0A=-1  $ 

91 

SAVE 

MPCT1, SINGLE, OMIT, W09ET, REACT, MPCF2,NSKIP, REPEAT, NOL, NO A  $ 

92 

OPDB 

EST  ,  OYNAMICS/RAD1  ,R!A02  yCONV,  HTGEN.FLUX/V  ,N,  NORAD  1/V  ,N  .N0RA02/  V, 
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RIGID  FORMAT  14 

N  A  STRAW  SOURCE  PROGRAM  COMPILATION 

ONAP-OMAP  INSTRUCTION 

NO. 

N,NOCON/V,N,NOGEN/V:,M, NOFLUX  $ 

93 

SAVE 

NORAOlt NOR A  02, NOCON, NOGEN, NOFLUX  $ 

94 

PURGE 

R02'T  AB,HTTAIB,FLXTABi/IMl/HTGE  N/NOGEN/CONVPT ,  CONT  AB,  KGGCON,  LOCON, 
KGSSUM/NOCON/  S 

95 

CHKPNT 

RADI  , RAD2 ,OONV, HTG0N.FL UX  t 

96 

LABEL 

LABL91  S 

97 

COND 

LABL73.NOCON  ’* 

98 

CONO 

LABL  81, IFIRST  S 

99 

CONV1 

GPCT  ,CONV/CQNVPT,CONTA&  S 

100 

CHKPNT 

CONVPT.CO'NT  AB  S 

101 

LABEL 

LABL  81  $ 

102 

C0NV2 

GPCT  , CONVPT.OIT, CONT AB.CONV/ KGGCON, LDCON/V,  N, T I ME/ V , N ,L USET/ V , 
N.STEP  S 

10  3 

ADO 

KGS,  KGGCON/ KGGSUM/C1,N,  11. 0,0.0) /C,N, 11.0,0.0)  $ 

10  4 

EQUIV 

KGSXY.KGG/IP1S 

105 

EQUI  V 

KGSS  UM, KGG/ 1  Ml  S 

106 

LABEL 

L AB L  73  S 

10  7 

PURBE 

GM.'G  MD/MPCF 1  /  GO ,  KOOIB,  LOO  ,UOO,MOOB,MOAB,GOD/OMIT/KFS,PST,QP/ 
SINGLES 

10  8 

EQUIV 

KGS , KNN/WPCF 1/MGG,MNN/MPCF1/  BGG , BNN/MPCF1/ K4GG, K4NN/MPCF1S 

109 

CHKPNT 

GM,RG,GO, KOOB, LOO, UOU, NOOB, MOAB , KFS,QP, USET, GMD, GOD, PST,  KNN , 
MNN , BNN.K4NN  S 

110 

CONO 

LB14,6ENEL  S 

111 

CONO 

LBL4.N0SIMP  S 

112 

GPSP 

GPL,  GPST,  USET,  SIL/tTGPST  S 

113 

OFP 

OGPST,,,,|,//V,N , CARIONO  S 

114 

SAVE 

CARO  NO  S 

3.15-5  (7/1/72-NSRDC) 


RIGID  FORMAT  14 

N  A  S  T  R  A  N  SOURCE  PROGRAM  COMPILATION 
OMAP-OMAP  INSTRUCTION 

NO. 

115 

label 

LSI  4  $ 

116 

CONO 

LBL2.MPCF1  S 

117 

COND 

LBL  220» IFIRST  S 

116 

MCE1 

USET.RG/GM  $ 

119 

CHKPNT 

GM  $ 

120 

LABEL 

LBL  2 20  S 

121 

MCE2 

USET  ,GM,KGG,MGG,BGG,K4QG/KNN,MNN,8NN,K4NN  0 

122 

CHKPNT 

KNM  , MNN»BNN, K4NN  S 

123 

LABEL 

LBL  2  $ 

124 

EQUI  V 

KNN  ,  KFF/SINGLE/MNN.IMFF/S INGL  E/BNN*  BFF/SI NGL E/K4NN,  K4FF/SINGLES 

125 

CHKPNT 

KFF, MFF.8FF.K4FF  $ 

126 

COND 

L8L3, SINGLE  $ 

127 

SCE1 

USET  ,KNN,MNN,8NN,K4>NM/KFF,KFS,KSS,MFF,BFF,K4FF  X 

128 

CHKPNT 

KSS  S 

129 

CHKPNT 

KFS  ,  KFF,MFF ,  8FF ,  K4FIF  S 

130 

LABEL 

LBL  3  $ 

131 

EQUIV 

KFF  ,  KAA/OMIT/  MFF  ,M  AAl/CJMIT/BFF,  B  A  A/OMIT /K4FF ,  K4  A  A/  OMI  T$ 

132 

CHKPNT 

KAA, MAA,BAA,K4AA  $ 

133 

CONO 

LBL5.0MIT  $ 

134 

SMP1 

USET  ,KFF,MFF,BFF,K4<FF/GO,KAA,KOOB,LOO,UOO,NAA,MOOB,MOA8,BAA, 

K4A  A  S 

135 

CHKPNT 

G0*KAA,MAA,BAA,K4AA  *• 

136 

LABET 

LBL5  $ 

137 

EQUI  V 

KAA.,  K00/IM1  S 

138 

CONO 

LBL  2  30, N0K4GG  $ 

139 

ADD 

BA A , K4A A/BDD/C, N, (1.0,0.0)/C,N,  (1.0, 0.0)  S 

3. 15-6  (7/1/72-NSRDC) 


• 

RIGID  FORMAT  14 

DMAP 

NO. 

N  A  S  T  R  A  N  SOURCE  WWOGRAM  COMPILATION 
-OMAP  INSTRUCTION 

140 

JUMP 

LBL231  $ 

141 

LABEL 

LBL230  $ 

142 

EOUI V 

BAA, BOO /I Ml  S 

143 

LABEL 

L8L231  S 

144 

CHKPNT 

KOO , 800  t 

145 

COND 

LBL 108, IPIRST  S 

1V6 

OPD 

OYNAMICS,GPL,SIL,US£TVGPLD,SILD,USETO,,OLT,  ,  ,NLFT,TRL,,EQDYN/V, 

N, LUSET /V ,N, LUSETO^V,  N**NOTFL/V, N, N00LT/V,N, NOPSOL/ V ,N  NOFRL/V, 
N,N0NLFT/V,N,N0TRL/'V,NVN0EED/C,N,123/V,N,N0UE  s 

14  7 

SAVE 

LUS€TO,NOOLT,NONLFr,NOTRL,NOUE  S 

148 

PURGE 

PNLT)/NONLFTS 

149 

EQUIV 

GO jG  OD/NOUEVGM, GMO^NOUE  S 

ISO 

CONO 

LBL221.IM1  S 

1S1 

TABPT 

GO, XAA,BAA,CASEXZ,KGQ  //  $ 

• 

152 

LABEL 

LBL  221  S 

15  3 

CHKPNT 

USETO.EQOYM.OLT.TRL!, GOO, GMD,NLFT,PNLD, SILO,  GPLD  t 

154 

CONO 

ERR'ORI.NOTRLJ 

155 

PARAM 

//C,N,AOO/V,N,  NE  VERJ/O,  N,  1/C  ,  N ,  0  $ 

156 

PARAM 

//C  ,  N,MPY/V,N,REPEATnVC,N,l/C,N,-l  $ 

157 

JUMP 

LBL 13  S 

158 

LABEL 

LBL  1  3$ 

169 

PURGE 

PNLO,OUDV1,OPNL1,OUOV2,OPNL2,XYPLTTA,OPP1,OQP1,OUPV1,OES1, 
OEF1,OPP2,OQP2, 0UPV2 ,0ES2,0EF2,PL0TX2,XYPLTT/NEVER  S 

160 

EQUIV 

CAS£CC,CASEXZ/IM1S 

161 

CASE 

CASECC,  /CASE  XX/C  , N,IT(TA N/  V ,N,  REP E  ATT /  V,  N,  NOL  OOP  S 

162 

SAVE 

REPEATT, NOLOOP  S 

15  3 

CHKPNT 

CASEXX  s 

• 

l _ 

3.  15-7  (  7/1/72-NSRDC) 

RIGID  FORMAT  14 

N  A  S  T  R  A  *  SOURCE  WHO  G  R  AM  COMPILATION 
OMAP-DMAP  INSTRUCTION 


164  OPOC  USET  D,KFS, YS  jGM/REOUCE.'KFSYS/V*  N»MPCFi/V»N»  SINGLE/ V,N  .LUSETD/V, 

N,r$ 

165  CHKPNT  REDUCE  S 

166  PURGE  MATTAB,NMATAB,£CPT1,ECPT2/MAT/DLT/N0DLT/NLFT/N0NLFT/GM/MPCF1/ 

YSiKFSYS/Y  S 

167  DPOD  RAI1,RA02,HTGEN,CASEOC,<TRL,FLUX/RD1TAB,RD2TAB,HTTAB,TICDAT, 

TSrEP0,FLXTAB/V,N,SiCA».AR/V,N,N0RA01/V,N,N0RAD2/V,N,N0GEN/V,N, 
LUD/V,N,LLOAON/V,N,1NLFTP/V,  N, NOFLUX  $ 

16S  SAVE  SCALAR, NORA01,NORAO!2, NOGEN, LUD,LLOAON,NLFTP,  NOFLUX  S 

169  CHKPNT  RD1T AB,  RO  2T  AB,  HTTAB.,TICOAT,  TSTEPD  ,FLXTAB  S 

170  EQUIV  CASEXZ.CASECC/IPl  f 

171  LABEL  LBL108  S 

172  TROTEMP  OLT«,  NLFT,  OIT  ,  KOO , BUD, REDUCE ,  LOCON, EST  ,ECPT,  MATTA  B,  RD1 TA  B, 

TICD  AT,  TSiTEPO , KFS ,GH,USETD ,  YS*KFSYS/UDVT  ,  PPT  , PNL 0  .EST  Y  ,ECPT  Y, 
ECP7  1, ECPT2,  NMAT AB.JTBMD AT/V,  N,MAT/V,N, SC ALAR/V,N,TIME/V,N, STEP/ 
V,l»,  NOCON/\A,N,NOGEN/V,N,NORADl/V,N,NORAD2/V,N,LUSETD/V,N, 
SilfGLE/ V,  N,  MPCF1/V,JN»‘L17D/ V,  N,KCON/V»N,KCAP/V,N,IECPTY/V,N* 
IREPET/V,N,OMEGA/V»Y,£PSCON=0.001/V, N, IGROUP/V ,N,ISTEP/V,N, 

NSTE  P/V,N,WOUTA/V,*T,LLOADN/V,N,NLFTP/V,N,  NOFLUX/ V,Y,BT$ 

17  3  SAVE  MAT',  SCALAR,  TIME,  STEV,  NOiCON,  NOGEN,  NORAD1,  NOR  AD2,LUSETD, SINGLE, 

MPCF1,LUD,KCON,KCAP!,IEOPTY,IREPET,OMEGA,EPSCON,IGROUP,ISTEP, 
nstep,nouta,lloadn,inuftp,noflux,bt  $ 

174  CHKPNT  UDVT  ,  PPT,  PNL  0  ,ESTY  ,fECP?Y  ,  ECPTl,  ECPT2  ,  NMATAB  ,  TEMO  AT  $ 

175  EQUIV  EST  Y »EST/ I Wi/ECPTY, ECPT/ I Ml  $ 

176  EQUIV  NMAT  AB,MATTAB/IM1/TEMDAT ,TIC0AT/IM1  $ 

177  EQUIV  EST  , EST Y/ IP1 /ECPT , SEPT*/ IPi  S 

178  EQUIV  MATT  A8,  NM  AT  A  B/IPl/T'ICDAT  ,TEM0AT/IP1  S 

179  CONO  LBL  8 10, IM1  $ 

180  TABPT  ESTY,ECPTY,MATTAB,THODAT,//  S 

181  TABPT  NMAT  AB, TEMO AT ,KGG,CASEXZ,//  S 

182  LABEL  LBL810  $ 


3. 15-8  ( 7/1/72-NSRDC) 


RIGID  FORMAT  14 

NASTRAN  SOURCE  PROGRAM  COMPILATION 
OMAP-DMAP  INSTRUCTION 


NO. 

183 

COND 

LBL815« IFIRST  $ 

164 

MOOA 

/, * « /C,N,0.0/C,N,Q.U/C,N,0.0/C,N,0.Q/C,N,0. 
C,N,0/C,N,0/C,N,0.0/V,N, IFI R/V, N , IFIRST  t 

0/C,N,0/C,N,0/C,N,0/ 

185 

SAVE 

IFIR, IFIRST  $ 

166 

LABEL 

LBL  8  15  S 

187 

CONO 

LABL  74* IREPET  t 

188 

CONO 

LABL  95* IECPT  Y  $ 

189 

REPT 

LAB L  83. 50  0  S 

190 

JUMP 

ERRI0R4  S 

191 

LABEL 

LABL  95  S 

192 

COND 

LBL310.MAT  S 

193 

PARAM 

//C  ,  N,ANO/V,N,NBOTU/V:,N,  KCON/V,N,KCAP  $ 

194 

COND 

LBL  3 10  *  NBOTH  $ 

195 

JUMP 

LBL311  S 

196 

LABEL 

LBL'310  t 

197 

EQUIV 

KGS  X  Y  t  KGGX/ I  Ml  t 

198 

JUMP 

LABL  97  $ 

199 

LABEL 

LBL  311  * 

20  0 

COND 

LA8L96.KC0N  $ 

201 

S  MAI 

CSrM,MPT,ECPTl,GPCr,Dn*/KGGl,K4GGl,GPSTl/V, 
N0K4GG  t 

N, NOGENL /V ,N • 

202 

SMA1 

CSTM,MPT,ECPT2,GPCr.,DIT/KGG2,K4GG2,GPST2/V, 
N0K4GG  $ 

N, NOGENL/V ,N , 

20  3 

ADD 

KGS2  ,  KGG1/KGG3/C  »N»'(1.  0,0.0)  /C»N»  (-1.0, 0.0) 

$ 

204 

ADO 

KGSXY,KGG3/KGG4/C,W»*i.:0,0.0>/C,N,  (1.0, 0.0) 

s 

20  5 

EQUIV 

KGS4 , KGGX/I Ml  S 

206 

COND 

LABL 96, N0K4GG  S 

3. 15-9  (7/1/72-NSRDC) 


RIGID  FORMAT  14 

NASTRAW  SOURCE  PiRIOGRAM  COMPILATION 
DM AP-OMAP  INSTRUCT  ION 


NO. 

207 

ADD 

K4GG  2,K4GG1/  K4GG3/CI*  N,  (1.0,0.01  /C, N, (-1 . 0 ,0 . 0 )  $ 

208 

ADD 

K4SG,K4GG3/K4GG4/C,IN,'  (1.0.0. 0>/C,N,  (1.0 , 0.0)  S 

209 

EQUIV 

K4G G4.K4GG/ 1 Ml  S 

210 

LABEL 

LABL96  S 

211 

CON3 

LABL97, KCAP  S 

212 

SMA2 

CSTM,MPT,ECPT1,GPCT*,0IT/,BGG1/V, Y.WTMASS/V, N,NOMGG/V, N, NOBGG/V, 
Y.COUPMASS  S 

213 

SAVE 

NOMGGtNOBGG  $ 

214 

SMA2 

CSTM,MPT,'ECPT2,GPCT<,01T/,BCG2/V,  Y, WTMASS/V, N ,NOMGG/ V, N, NOBGG/ V, 
Y,  C'OUPMASiS  $ 

215 

SAVE 

NOHGG, NOBGG  S 

216 

ADD 

BGK2 ,  8GG1/BGG3/C *N,M1.0»0.0>/C,N, (-1.0,0.01  $ 

217 

ADO 

BGG , BGG3/BGG4/C, N, Cl. 0,0 .0) /C,N, (1.0, 0.0)  $ 

218 

EQUIV 

BGC4  ,BGG/IM1  S 

219 

LABEL 

LABL97  S 

220 

REPT 

LABL  98, 50 0  I 

221 

TABPT 

DLT  ,  NLFT»iB0D»RADl,RIA02//  S 

222 

TABPT 

HTSEN, REDUCE, LDCON.IRCTITA  B.RD2TAB//  $ 

223 

TABPT 

HTT  A  B.TICOAT , TSTEPD.FLUX *ELXT  AB//$ 

224 

TABPT 

COMV  PT ,CONT  AB,CONV,XGGV//$ 

225 

JUMP 

ERRI0R4  S 

226 

LABEL 

LABL74  $ 

227 

EQUIV 

PPT.PDT/NOSET  S 

228 

CHKPNT 

UOVT,PDT,PST,PPT,P(fLO  S 

22  9 

VDR 

CAS£XX,EQDYN  ,USETD,'UOVT  ,  PPT  ,  XYCD  B.PNLD/OUOV  1 ,  OPNLi/  C,N, 

TRA NRESP/C, N *  DIRECTVC.N, 0/V,N,NOO/V, N, NO P/C , N, 0  $ 

230 

SAVE 

NOD, NOP  $ 

3. 15-10  (  7/1/72-NSRDC) 


RIGID  FORMAT  14 

NAS  TRAM  SOURCE  PROGRAM  COMPILATION 
ONAP-DMAP  INSTRUCTION 


NO. 

231 

CHKPNT 

OUOVltOPNLl  S 

232 

COND 

LBL15.N0D  S 

233 

SDRS 

OUOVltOPNLl,  ,.,/OUD-V2,OPNL2,  ...  $ 

234 

OFP 

OUO  V 2t OPNL2  ,  •  •  ,//V»IN,CARDNO  t 

235 

SAVE 

CA'RDNO  t 

236 

CHKPNT 

OPML 2,OUOV2  $ 

237 

XYTRAN 

XYCTO  B,  OUOV2,  OPNL  2,,  l.^XYPLTT  A/C,  N  ,  TRAN/C,  N,DSET/V,N,  PFILE/  V. 
CARO  NO  t 

238 

SAVE 

PFILE, CARONO  t 

239 

XYPLOT 

XYPLTTA//  S 

240 

LABEL 

LBL15  S 

241 

PARAM 

//C\  N,  AND/V,  N, PJUMPVY,  M,  NOP/ V,  N,  JUMP  PLOT  $ 

242 

COND 

LBL 1 8 ,P JUMP  t 

243 

EQUIV 

UOVT,UPV/NOA  % 

244 

COND 

LBL17.N0A  t 

245 

SDR1 

USET 0 ,  ,UO VT  ,  ,  YS  ,  GOD!,  ©My  ,KFS  , KSS  , /UPV,  ,QP/C ,  N,  1/C,  N,  DYNAMICS 

246 

LABEL 

LBL17J 

247 

CHKPNT 

UPV, OP  $ 

248 

S0R2 

CASE XX.CSTM, MPT ,DIT<,E:QDYN,SILD»  ,  *BGPDT,PPT ,  QP,UPV,EST,XYCDB/ 
OPPtl  ,OQPl,OUPVl,OESil,OEFl,PUGV/C,N, THERMAL  S 

249 

50R3 

OPPl,OOPl,OUPVl,OESa,OEFl,/  OPP2,OQP2,OUPV2,OES2,OEF2,  $ 

250 

CHKPNT 

OPP12  ,OQP2,OUPV2  ,  OES2,OEF 2  S 

251 

OFP 

OPP2  , OQP2  ,OUPV2 , 0EFI2.0ES2  ,//V,N,CARDNO  $ 

252 

SAVE 

CARONO  S 

253 

COND 

P2,  JUMP  PLOT  S 

254 

PLOT 

PLTP AR, GPSETS .ELSETStCAS  EXX , BGPDT ,EQEXIN,SIL,, PUGV/ PL OT X2/  V 

NSI L/V, N,  LUSET/V  ,N»(JUMPPLOT /  V»N , PLTFLG/V ,N, PFILE  S 


3. 15-11  ( 7/1/72 -NSRDC) 


RIGID  FORMAT  14 

N  A  S  T  R  A  N  SOURCE  PROGRAM  COMPILATION 
OMAP-OMAP  INSTRUCTION 


NO. 

255 

SAVE 

PFILE  S 

256 

PRTMSG 

PLOT  X2//  $ 

257 

LABEL 

P2  S 

258 

XYTRAN 

XYCOB.OPP2,OQP2,OUPIV2,OES2,OEF2/XYPLTT/C,N,TRAN/C.N,PSET/  V.N, 
PFILE/V.N.CARDNO  $ 

259 

SAVE 

PFU.E.CARONO  S 

260 

XYPLOT 

XYPLTT//  $ 

26-1 

LABEL 

LBL18  $ 

262 

PURGE 

SK6GX,SGPST,SMGG,SOGPMQ,SKGG,SRG,SYS,SUSET,SOGPST/IFIR  S 

263 

PURGE 

SGH,SKNN,S*FF,SKFS,ISKS9,SGO,SKAA,SKOOB,SLOO.SUOO,SKLL/IFIR  $ 

264 

PURGE 

SKLR,SKRR,ST.LL.SULU.SOH,SPG,SCIR,SPO,SPS.SPL/IFIR  $ 

265 

PURGE 

SUL  V ,  SUOOV , SRUL V .  SRfUOV , SUGV » SPGG  .SQG.0PG1.0QG1.0UGV1/IFIR  t 

266 

PURGE 

0ES61,0EFF1,PUGV1,0€SAVG,PLTX2,CASEXY,GE0M3X,ESTZ,ECPTZ/IFIR  S 

267 

PURGE 

GPTTX.GEOM3/IFIR  S 

268 

JUMP 

LABL99  S 

269 

LABEL 

LABL99  t 

270 

RETEMP 

CAS£XZ,GEON3,UPV,E»OYN,EQEXIN,SILD,EST,ECPT/CASEXY,GEOM3X,ESTZ, 
ECPtTZ/V,N,NSTEP/V.W,»GAIN/V»N,NREC/V,N,  NOSUB/V  ,  N  ,  IC  ASE  S 

271 

SAVE 

NSTEP, AGAIN*  NREC.NC7SU6.  I CASE  S 

272 

CONO 

FINIS. NOSUB  $ 

273 

CHKPNT 

ESTZ.ECPTZ.CASEXY  $ 

274 

EQUIV 

EST<Z  .EST/IM1/ECPTZ.ECPT/IM1/CASEXY .CASECC/I  Ml  $ 

275 

EQUIV 

EST  .  ESTZ/IP1/ECPT  .ECPTZp/  IPi/CASECC.C  ASEXY/IP1  % 

276 

CONO 

LBL111.IM1  | 

277 

TABPT 

ECPTZ,ESTZ,CASEXY,,l  //  S 

278 

LABEL 

LBL111  r 

3.15-12  ( 7/1/72 -NSRDC) 


RIGID  FORMAT  14 

NASTRAN  SOURCE  PI  R  8  G  R  AH  COMPILATION 
ONAP-DNAP  INSTRUCTION 


NO. 

279 

GP3 

GEBM3X,EQEXIN,GE0M2>,/  .GPTTX/C  ,N,123/V,N,  N06RA  V/C  ,N  t123  S 

288 

SAVE 

NOSRAV  S 

281 

EQUIV 

GPTT  X.GPTT/IM1  S 

282 

CHKPNT 

GPFTS 

283 

CONO 

LBL201. NOSIMP  * 

284 

CONO 

L8L909, MOMATt 

285 

PARAN 

//C*.  N,SU8/V,N,IDIFF1/¥,N,  ICASE/V.  N.IFOURS 

286 

COND 

LBL909* IOIFFS 

28  7 

JUMP 

LBL'201S 

28  8 

LABEL 

LBL909S 

289 

SMA1 

CST»M  ,  MPT*  ECPT «  GPCT'.IOST/SKGGX  »  .SGPST/V,  N»  NOGENL/V  ,N, N0K4GG  S 

290 

CHKPNT 

SGPST.SKGGX  S 

291 

CONO 

LBL  201*  SKPMGG  8 

292 

SMA2 

CSTM  .MPT, ECPT .GPCT ,OIT/'SMGG»  /V,  Y.NTMASS/V.N,  NOMGG/ V ,N, NOB GG/V, 
Y, C  OUPMASS  S 

293 

SAVE 

NOMGG  t 

294 

CHKPNT 

SMOG  S 

295 

COND 

LBL  201,  GRIOPNT  S 

296 

CONO 

ERROR7, NOMGG  S 

297 

GPHG 

BGP OT, CSTM, EQEXIN.SlMGG/SOGPHG/V,  Y,GRDPNT*-1/V,Y,  WT MASS  S 

298 

OFP 

S08PHG,  ,*  ,,//V,N,CA.RDNO  S 

299 

SAVE 

CARO  NO  $ 

300 

LABEL 

LBL 201  S 

301 

EQUIV 

SK8GX.SKGG/N0GENL  S 

30  2 

CHKPNT 

SKGG  $ 

30  3 

COND 

LBL  211, NOGENL  S 

3.15-13  ( 7/1/72-NSRDC) 


RIGID  FORMAT  14 

N  A  S  T  R  A  M  SOURCE  PI  R  0  GR  AN  COMPILAT  I  ON 

ON AP-OHAP'  INSTRUCTION 

NO. 

38% 

SNA3 

GEI',  SKGGX/SKGG/V,N,1_USET/V,N,NOGENL/V,N,NOSIMP  $ 

309 

CHKPNT 

SKSG  S 

306 

LABEL 

LBL  2 11  S 

307 

PARAM 

//c,  n,mpy/f,n,nnskiip/cvn,o/c,n,o$ 

308 

GP4 

CASECC»GE0H4 , EQEXIW,SIL , GPDT/SRG , SYS. SUSET  ,/V»N»LUSET/V»N, 
MPCF1/V,N,NPCF2/V,N,SINGLE/V,N,0MIT/V,N,REACT/V,N,NNSKIP/V,N 
REPE  AT/V»iNf‘NOSET/V»lN»NOiL/V»  N,NOA  $ 

30  9 

SAVE 

MPCF 1,  MPCF2,  SINGLE  ,IONIT,  REACT, NNSKIP,  REPEAT,  NOSET.NOL.NO  A  S 

310 

CONO 

ERRORS,  NOL  S 

311 

PARAH 

//C, N , A  ND/V , N , NO  SR i^V  ,  N  ,  S INGL  E/  V  ,  N,  RE  ACT  S 

312 

PURGE 

SKRR,SKL*,SQR,SDM/RfE*CT/SGM/MPCFl/SGO,SKOOB,SLOO,SUOO,SPO, 
SUOOV,SRUOV/  OMIT/SPlS,  SKF  S»  SK  SS/S  INGL  E/SQG/N  OSR  $ 

313 

EQUIV 

SKKG.SKNN/HPCFl  S 

314 

CHKPNT 

SKRIR,  SKLR,  SQR,  SD  H, S.GH,  SGO,  SKOOB,  SL  00,  SUOO,  SPO.SU  00  V,SQG,SPS, 
SKFS  .SKSS1,  SUSET  .SRG.SYS,  SRUOV.SKNN  S 

315 

COND 

LBL  204, GENEL  $ 

316 

GPSP 

GPL, SGPST, SUSET, SIU/SOGPST  S 

317 

OFP 

SOtPST,,,,  ,//V,N,CAiRONO  $ 

318 

SAVE 

CARO  NO  S 

319 

LABEL 

LBL 204  S 

320 

COND 

LBL  2 13, MPCF2  1 

321 

MCE! 

SUSCT  ,SRG/S>GM  $ 

322 

CHKPNT 

SGR  $ 

323 

HCE2 

SUSET ,S6M,S,KGG,  ,  ,/SKNN,  ,  ,  $ 

324 

CHKPNT 

SKNN  S 

325 

LABFL 

LBL  2 13  $ 

326 

EQUIV 

skbn,skff/single  s 

3.15-14  (  7/1/72-NSRDC) 


RIGID  FORMAT  14 

nastram  source  program  compilation 

OMAP-DMAP  INS  TRUCd ION 


NO. 

327 

CHKPNT 

SKFF  $ 

328 

CONO 

L8L215* SINGLE  t 

329 

SCE1 

SQSET  tSKNN.  «  «/SKFF|lSKFS.SKSS  ...  S 

330 

CHKPNT 

SKFS.SKSS.SfKFF  $ 

331 

LABEL 

LBL215  S 

332 

EQUIV 

SKFF  »SKAA/OMIT  $ 

333 

CHKPNT 

SKAIA  $ 

334 

CONO 

LBL216.0MIT  $ 

335 

SMP1 

SUS£T,SKFF,  ,  ,/SGO,S.KAA«SKOOB,SLOO,SUOO,,  ,,,  $ 

336 

CHKPNT 

SGO.SKAAfSKOOBf  SL0O.5U00  $ 

337 

LABEL 

LBL216  S 

338 

EQUIV 

SKAiA  tSKLL/REACT  S 

339 

CHKPNT 

SKLL  S 

340 

CONO 

LB216, REACT  S 

341 

RBMG1 

SUS£T,SKAA,/SKLL.SKLRF»SKRR,»  ,  S 

342 

CHKPNT 

SKLL.SKLR.SKRR  S 

343 

LABEL 

LB216  S 

344 

RBHS2 

SKLL/SLLL.SULL  S 

345 

CHKPNT 

SULL.SLLL  $ 

346 

CONO 

LBL217, REACT  $ 

347 

R9MG3 

SLLL,SULL,SKLR,SKRRI/SDM  $ 

348 

CHKPNT 

SDN  $ 

349 

LABEL 

LBL217  t 

350 

SSG1 

SLT'»  BGPOT »CSTM»  SILflEST <MPT»  GPTT , EOT , SMGG .CASECC , DIT/SPG/V ,N 
LUS£T/V»N»NNSKIP  S 

351 

CHKPNT 

SPG  $ 
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RIGID  FORMAT  14 

NASTRAM  SOURCE  PROGRAM  COMPILATION 
OMAP-OMAP  INSTRUCTION 


NO. 

352 

EQUIV 

SPG.SPL/NOSET  $ 

353 

CHKPNT 

SPL  $ 

354 

COND 

LBL210. NOSET  $ 

355 

SSG2 

SUSET  fSGM.SYS.  SKFS.ISGJO  ^  SDH»  SPG/SQR*SPO»SPS»  SPL  t 

356 

CHKPNT 

SQR» SPO.SPS*  SPL  t 

367 

LABEL 

LBL210  $ 

358 

SSG3 

SLLL  tSULL.SKLL * SPL.iSLOO » SUQO « SKOOB.SPQ/SULV *  SUOO V, SRULV  »  SRUO V7 
V,N,OMIT/V,Y,IRES=-l/-V,N,NSKIP/V,N,EPSI  S 

359 

SAVE 

EPS!  t 

360 

CHKPNT 

SUL  V  .SUOOfV « SRULV  •  SRJUOV  S 

361 

CONO 

LBL9.IRES  $ 

362 

MATOPR 

GPL.USET,6IL,SRULV//0.N, L  S 

363 

MATGPR 

GPL.USET,SIL,SRUOV/!/C,N,0  S 

364 

LABEL 

LBL9  S 

365 

SOR1 

SUSET «SPG « STJL  V  •  SUOOIVt'SYS  •  SGO  *SGN»  SPS , SKFS.SKSS . SQR/SUGV  .  SPGG  , 
SOt/ V,N»NNS KIP/C  «N. STATICS  S 

366 

CHKPNT 

SU6V.SPG6.SQG  S 

36.7 

SDR2 

CAS ECC.CSHM. MPT, OITI,eQCXIN»  SIL.  GPTT,  EDT  ,  BGPDT, SPGG, SQG.SUGV, 
EST'»  /OPG1 , OQG1 , OUGVl.OESSl* OEFF 1  »PUGV1/C*N,  STATICS  $ 

368 

OFP 

OUt  V 1  *OPGl  t  OQG1*  OEFIF1  tOE  SSI  *  //V  »  N*CARONO  S 

369 

SAVE 

CARONO  $ 

370 

STRSAVG 

EOEX I N» 0ES31/0ES AVG  t 

371 

OFP 

OESJ4VG*  .«  **//V*  N,  CAIRO  NO  $ 

372 

SAVE 

CARONO  t 

37'3 

CONO 

P22'»  JUMPPLOT  S 

374 

PLOT 

PLTiPAR.GPSETS*  ELSETSi C ASECC » BGPDT  tEQEXINtSIL  t PUG VI ,/PLTX2/V»N» 
NSIL /V»  N»  LUSET /V  »N»  JUMPPLOT /V»N  »  PLTFLG/V*N. PFILE  S 

3. 15-16  (7/1/72-NSRDC) 


RIGID  FORMAT  14 

N  A  S  T  R  A  N  SOURCE  Pi  R  •  G  R  A  M  COMPILATION 
DMAP-OMAP  INSTRUCTION 


NO. 

3/5 

SAVE 

PFILE  t 

3/6 

PRTMSG 

PLTIX2//  f 

3// 

LABEL 

P22  $ 

3/8 

CONO 

FINIS. AGAIN  $ 

3/9 

REPT 

LABL99. 110  S 

380 

FILE 

SK*GX*SAVE/SKGG  =  SAVE/»SKNN=SAVE/SKFF*SAVE/StCAA*SAVE/SKLL*»SAVE 

381 

JUMP 

ERRIOR5  S 

382 

LABEL 

ERROR 1$ 

363 

PRTPARM 

//C,  N,-l/C,N,OIRTRD:$ 

384 

LABEL 

ERROR4  t 

385 

PRTPARM 

//C, N,-4/C,N,OIRTRB  % 

386 

LABEL 

ERRORS  S 

38/ 

PRTPARM 

//C',N,-1/C»N,ST  AT  ICS  t 

388 

LABEL 

ERROR6  S 

389 

PRTPARM 

//C",N, -4/C, N, STATICS  S 

390 

LABEL 

ERROR7  * 

391 

PRTPARM 

//C.N, -3/C, N, STATICS  S 

392 

LABEL 

ERRORS  S 

393 

PRTPARM 

//C, N, -2/C, N, STATICS  S 

394 

LABEL 

FINIS! 

395 

END 

s 
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3.15.2  Description  of  DMAP  Operations  for  Thermal  Transient - 
Structural  Statics  Combinations 

DMAP  statements  1-260  are  associated  with  the  thermal  transient 
portion  of  the  analysis,  while  DMAP  statements  261-379  are 
associated  with  the  structural  statics  portion.  Statements  380-395 
deal  with  error  conditions.  Only  selected  DMAP  statements  or  sets 
of  statements  will  be  described.  Descriptions  of  other  DMAP 
statements  may  be  found  in  the  Description  of  Static  Analysis, 

Section  3.2.2  and  Direct  Transient  Response,  Section  3.10.2. 

4-6.  PARAM  sets  up  DMAP  parameters  for  later  use. 

7.  MODA  initializes  DMAP  parameters  IFIR  and  IFIRST. 

MO  DA  is  later  used  to  switch  the  values  of  these  parameters. 
42.  DPDA  converts  temperature-dependent  thermal  material 

tables  into  step  functions  based  on  TRANGE  bulk  data  cards 
and  creates  a  new  DIT  data  block. 

50.  INITEM  modifies  the  EST  and  ECPT  data  blocks  by  inserting 

the  initial  element  temperatures  as  derived  from  the 
initial  conditions  specified  on  TIC  bulk  data  cards. 

INITEM  also  creates  data  block  MATTAB  which  contains 
the  temperature  interval,  as  derived  from  TRANGE  cards,  in 
which  each  thermal  material  property  for  each  element  falls. 
92.  DPDB  computes  area,  surface,  and  volume  integrals 

required  by  convection,  radiation,  and  heat  generation. 
99,102.  CONVl,  CONV2  compute  contributions,  due  to  convection. 

164.  DPDC  creates  a  matrix  which,  when  premultiplies  a  load 

vector,  reduces  the  vector  to  take  account  of  single  point 
constraints  and  multi-point  constraints. 


3. 15-17a  (7/1/72 -NSR DC) 


187. 

188, 189. 


193,219. 


220. 

270. 


272. 

370. 

378. 

379. 


Go  to  statement  226  if  temperature  history  is  complete. 

Go  back  to  statement  60  to  recompute  all  element 
conductance  and/or  capacitance  matrices  if  more  than  half 
the  matrices  must  be  changed. 

If  less  than  half  the  element  matrices  must  be  changed, 
recompute  the  necessary  matrices,  compute  the  net  effect 
of  changing  these  elements,  and  add  the  net  results  into  the 
structure  conductance  and/or  capacitance  matrices,  as 
needed. 

Go  back  to  statement  77  to  reduce  the  matrices  and  continue. 
RETEMP  interrogates  case  control  to  determine  a  time 
slice  selection.  If  one  is  found,  RETEMP  creates  TEMP 
bulk  data  cards,  recreates  GEOM3,  and  sets  counters  for 
the  next  entry  to  RETEMP. 

If  no  time  slice  selections  exist,  go  to  statement  394. 
STRSAVG  computes  average  grid  point  stresses. 

If  no  more  time  slice  selections  exist,  as  determined  by 
RETEMP,  go  to  statement  394. 

Go  back  to  statement  269  to  process  the  next  time  slice 
selection. 
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3.15.4 


Case  Control  and  Parameters  for  Thermal  Transient  - 


Structural  Statics  Combination 


The  following  items  relate  to  subcase  definition  and  data 
selection  for  Thermal  Transient-Structural  Statics  Combination. 

1.  If  more  than  one  subcase  is  defined,  the  first  subcase 
is  assumed  to  refer  to  the  thermal  transient  portion  of 
the  analysis,  and  all  other  subcases  are  assumed  to 
refer  to  the  structural  statics  portion.  If  one  or  no 
subcases  are  defined,  the  case  control  data  are  assumed 
to  refer  to  the  thermal  transient  portion  only,  in  which 
case  no  structural  static  analysis  will  be  performed. 

2.  After  the  first  subcase,  the  subcase  number  will  take 
on  a  special  meaning.  If  the  second  or  later  subcase 
number  is  N,  then  the  N-lst  output  time  step,  as  in¬ 
dicated  on  the  TSTEP  bulk  data  card  selected,  will  be 
the  time  for  which  the  static  analysis  will  be  performed. 

3.  DLOAD  or  NONLINEAR  must  be  used  in  the  first  subcase 

to  define  a  time-dependent  loading  condition  in  addition 
to  any  specified  thermal  loading  conditions. 

4.  TSTEP  must  be  used  to  select  the  time-step  intervals  to 

be  used  for  integration  and  for  output  for  the  first 
subcase.  The  time  steps  output  are  those  indicated  by 
the  skip  factors  on  the  selected  TSTEP  bulk  data  card 
and  the  last  time  for  each  specified  interval  (if  that  time 
was  not  already  indicated  by  the  skip  factors) . 
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5.  A  s.eparate  subcase  must  be  defined  for  each  static  analysis 
desired,  with  the  subcase  number  taking  on  the  meaning 
described  in  item  2. 

6.  A  static  loading  condition  will  automatically  be  defined  for 
each  static  subcase  by  NASTRAN.  This  condition  will  consist 
of  the  loads  induced  by  the  temperature  at  that  time.  (This 

is  equivalent  to  a  TEMPERATURE  (LOAD)  selection.)  Additional 
loads  may  be  defined  with  a  LOAD  selection  or  through  grid-point 
displacements  on  SPC  cards. 

7.  An  SPC  set  must  be  selected  for  each  subcase,  unless  the  model 
is  a  properly  supported  free  body,  or  unless  all  constraints 
are  specified  in  GRID  cards  or  on  Scalar  Connection  cards. 

8.  If  nonzero  initial  conditions  are  desired,  IC  must  be  used  to 
select  a  TIC  card  set  in  the  Bulk  Data  Deck.  Unlike  in  other 
rigid  formats,  TIC  cards  should  be  used  for  all  grid  points, 
not  just  those  on  the  analysis  set. 

The  following  printed  output,  sorted  by  point  number  or  element 
number  (S0RT2) ,  is  available  at  the  selected  output  time  steps  of 
the  thermal  transient  portion  of  the  rigid  formats: 

1.  Temperatures  (DISPLACEMENT)  and  rates  of  change  of  temperature 
with  respect  to  time  (VELOCITY)  for  a  list  of  PHYSICAL 

points  (grid  points)  or  SOLUTION  points  (independent  degrees-of- 
freedom) . 

2.  Nonzero  components  of  the  applied  load  vector  and  single  point 
forces  of  constraint  for  a  list  of  PHYSICAL  points. 

3.  Nonlinear  force  vector  for  a  list  of  SOLUTION  points. 
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The  following  plotter  output  is  available  for  the  thermal 
transient  portion  of  the  analysis: 

1.  Undeformed  plot  of  the  structural  model. 


2.  X-Y  plot  of  the  temperature  (component  1) 

or  rate  of  change  of  temperature  with  respect  to  time  of 
a  PHYSICAL  or  SOLUTION  point. 

3.  X-Y  plot  of  component  1  of  the  applied  load  vector  or 
nonlinear  force  vector- 

The  following  output  may  be  requested  for  the  static  analysis 
portions  of  the  rigid  format- 

1.  Displacements  and  nonzero  components  of  static  loads  and 
single  point  forces  of  constraint  at  selected  grid  points 
or  scalar  points. 

2.  Forces  and  stresses  in  selected  elements. 

3.  Undeformed  and  deformed  plots  of  the  structural  model. 
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The  following  parameters  may  be  used: 


1.  GRDPNT  -  optional  -  A  positive  integer  value  of  this 
parameter  will  cause  the  Grid  Point  Weight  Generator 
to  be  executed  and  the  resulting  weight  and  balance 
information  to  be  printed. 

2.  WTMASS  -  optional  -  The  terms  of  the  mass  matrix  are  mul¬ 
tiplied  by  the  real  value  of  this  parameter  when  they 
are  generated  in  SMA2. 

3.  COUPMASS —  optional  -  A  positive  integer  value  of  this 
parameter  will  cause  the  generation  of  coupled  mass  and 
capacitance  matrices  rather  than  lumped  mass  and  capacitance 
matrices  for  the  isoparametric  elements. 

4.  IRES  -  optional  -  A  positive  integer  value  of  this  para¬ 
meter  will  cause  the  printing  of  the  residual  vectors 
following  the  execution  of  SSG3  in  the  static  analysis 
portion  of  the  rigid  format. 

5 .  BT  -  required  if  radiation  is  being  used  and 

temperatures  are  not  specified  in  °R  or  The  usual 

values  for  BT  are  459.69  and  273.16  for  specifications  in 

°F  and  °C ,  respectively. 

6.  METCON  -  Optional  for  Rigid  Format  14.  A  value  of  1  will 
cause  the  program  to  make  convection  contributions  only  to 
the  load  vector.  For  full  details  see  Remark  4,  page  2.4-42c. 
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Some  other  points  to  be  noted  are: 


1.  The  only  degree-of-freedom  in  the  thermal  transient  portion  is  temperature/ 
which  is  assumed  to  be  component  1.  NASTRAN  automatically  contrains 
components  2-6.  Therefore,  SPC,TIC,  etc,  specifications  should 

contain  only  component  1. 

2.  Permanent  constraints  placed  on  a  GRDSET  card  or  on  GRID  cards 
will  be  used  for  the  statics  portion  only  and  are  ignored  for  the 
thermal  portion. 

3.  No  extra  points  or  direct  matrix  input  are  allowed. 

4.  TEMP,,  TEMPD,  TEMPP1 ,  TEMPP2 ,  TEMPP3 ,  and  TEMPRB  hulk  data  cards 
are  not  allowed. 

5.  All  OMIT  cards  will  be  ignored  in  the  thermal  portion. 

6.  All  CONVEC,  HFLUX,  HGEN,  RADIAT1,  and  RADIAT2  bulk  data  cards  present 
will  be  used  to  compute  thermal  loads. 

7.  Element  temperatures  will  always  be  used  to  compute  temperature- 
dependent  thermal  and  structural  material  properties.  (This  is 
equivalent  to  a  TEMPERATURE  (MATERIAL)  case  control  selection.) 

8.  In  the  statics  portion,  distributed  pressure  loads  specified  by 
the  PL0AD2  card  can  be  applied  only  to  elements  IS2D4,  SURF4, 

IS2D8,  and  SURF 8. 

9.  MATT  7,  MATT 8,  and  MATT 9  cards,  must  be  accompanied  by  corresponding 
TRANGE  cards.  See  the  discussion  of  the  TRANGE  bulk  data  cards. 

10.  On  restart,  no  rigid  format  switching  from  rigid  format  14  to  any 
other  rigid  format  or  from  any  other  rigid  format  to  rigid  format 
14  is  allowed. 
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11.  Consistency  of  units,  especially  for  temperature,  must  be  maintained. 
If  the  units  of  thermal  material  constants  are  per  °F  or  per  °R, 

all  temperatures  specified,  whether  through  SPC  cards,  TABLET1  cards, 
etc.,  must  be  in  °F  or  °R.  Similar  rules  apply  to  °C  and  °K.  The 
program  makes  no  conversions  for  differing  units.  Therefore,  the 
user  is  responsible  for  maintaining  consistency.  Input  and  output 
temperatures  will  be  in  the  same  units  for  the  statics  portion  of 
the  analysis,  so  care  must  be  taken  that  the  thermal  expansion 
coefficient  vector  (on  materials  cards)  is  given  in  the  correct 
units . 

12.  Time-dependent  temperatures  may  be  specified  at  a  grid  point  by 
connecting  a  scalar  element  from  the  grid  point  (degree-of-freedom  1) 
to  ground  with  a  very  large  "spring"  constant,  say  K^.  Then,  using 

a  TL0AD1  bulk  data  card,  specify  a  load  of 

P(t)  =  TK 

o 

where  T  is  the  desired  temperature  at  time  t.  If  is  large 
enough,  any  other  contributions  of  conductivity  or  load  at  that 
grid  point  will  essentially  be  neglected. 

13.  The  only  elements  presently  allowed  in  rigid  format  14  are  IS2D4, 
IS2D8 ,  IS3D8 ,  IS3D20,  SURF1,  SURF4,  SURF8,  and  BAR. 

14.  The  HEAT  parameter  option  must  not  be  used  with  rigid  format  14. 

15.  The  thermal  subcase  must  always  be  present,  even  on  restarts.  If,  on 
restart,  NASTRAN  determines  that  the  thermal  portion  was  completed 

on  the  previous  run  and  that  no  changes  have  been  made  to  the  thermal 
portion,  the  thermal  portion  will  not  be  re-executed. 
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598 


*** 


USER  FATAL  MESSAGE  998,  APPISO  CARD  REQUIRED 


If  isoparametric  elements  are  used,  the  APPISO  bulk  data 
card  must  be  present. 


999  ***  user  FATAL  MESSAGE  999,  NEITHER  TEMP  NOR  TEMPO  BULK  DATA 

CARDS  ARE  ALLOWED  IN  A  THERMAL  ANALYSIS  PROBLEM 

TEMP  and  TEMPD  bulk  data  cards  are  not  allowed  in  Rigid 
Format  14  (thermal  transient  analysis)  or  in  Rigid  Format  1 
(thermal  static  analysis) . 

999A  ***  USER  FATAL  MESSAGE  999A,  EPOINT  CARDS,  I.E.,  EXTRA  POINTS 

ARE  NOT  ALLOWED  IN  A  THERMAL  TRANSIENT  PROBLRM 

999B  ***  USER  FATAL  MESSAGE  999B,  TEMPP1,  TEMPP2 ,  TEMPP3 ,  AND  TEMPRB 

BULK  DATA  CARDS  ARE  NOT  ALLOWED  IN  RIGID  FORMAT  14 


1000  ***  USER  FATAL  MESSAGE  1000,  EITHER  THERMAL  TRANSIENT  OR 

THERMAL  STATIC  ANALYSIS  WAS  SELECTED  BUT  RIGID  FORMAT  14 
OR  1,  RESPECTIVELY,  WAS  NOT  SELECTED 

The  type  of  analysis  is  specified  on  the  APPISO  bulk  data 
card.  If  thermal  transient  analysis  is  selected,  Rigid 
Format  14  must  be  used.  If  thermal  static  analysis  is 
selected,  Rigid  Format  1  must  be  used. 

1000A  ***  USER  FATAL  MESSAGE  1000A,  HEAT  OPTION  AND  APPISO  BULK 

DATA  CARD  ARE  INCOMPATIBLE. 
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5151 


+  ** 


USER  FATAL  MESSAGE  5151,  MATT?  CARD  REFERENCES  UNDEFINED 
MAT 7  ****  CARD 


The  user  should  check  that  all  MATT7  cards  reference 
MAT7  cards  that  exist  in  the  Bulk  Data  Deck. 


5152  ***  USER  FATAL  MESSAGE  5152,  MATT8  CARD  REFERENCES  UNDEFINED 

MAT8  ****  CARD 


The  user  should  check  that  all  MATT8  cards  reference 
MAT8  cards  that  exist  in  the  Bulk  Data  Deck. 


5153  ***  USER  FATAL  MESSAGE  5153,  MATT9  CARD  REFERENCES  UNDEFINED 

MAT9  ****  CARD 

The  user  should  check  that  all  MATT9  cards  reference 
MAT9  cards  that  exist  in  the  Bulk  Data  Deck. 


5154  ***  USER  FATAL  MESSAGE  5154,  ELEMENT****  REQUIRES  AN  APPISO 

CARD 


An  APPISO  bulk  data  card  is  required  if  isoparametric 
elements  are  being  used. 


6.2-41  (7/1/72-NSRDC) 


5155 


*** 


USER  FATAL  MESSAGE  5155,  COORDINATE  SYSTEM  ****  IS 
NOT  A  RECTANGULAR  COORDINATE  SYSTEM 


The  material  coordinate  system  for  a  solid  isoparametric 
element  as  specified  on  a  CIS3D8  or  CIS3D20  bulk  data  card, 
must  be  defined  as  a  rectangular  coordinate  system. 


5156  ***  USER  FATAL  MESSAGE  5156,  ELEMENT  ****  PRODUCES  A 

SINGULAR  JACOBEAN 

Isoparametric  elements  require  that  the  determinant  of  a 
Jacobean  be  computed.  Check  the  geometry  for  the  element. 


5157  ***  USER  WARNING  MESSAGE "5157,  IN  MODULE  OESAVG  THE  ELEMENT 

TYPE  IS  ****  WHILE  THE  NUMBER  OF  WORDS  ON  THE  STRESS  FILE 
IS  **** 

STRESSES  WILL  BE  PRINTED  BUT  NOT  INCLUDED  IN  THE  AVERAGES. 

Logic  error  in  module  STRSAVG  or  averaging  of  stresses  of 
elements  with,  different  types  of  stresses  is  being 
attempted . 
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5158 


*  ** 


USER  WARNING  MESSAGE  5158,  IN  MODULE  OESAVG  TWO  DIFFERENT 
STRESS  COORDINATE  OUTPUT  SYSTEMS  HAVE  BEEN  SPECIFIED  FOR 
ELEMENTS  ****  AND  ****  BUT  THESE  ELEMENTS  HAVE  A  COMMON 
GRID  POINT. 

THE  STRESSES  FOR  THE  SECOND  ELEMENT  WILL  NOT  BE  INCLUDED 
IN  THE  AVERAGE 


The  stress  coordinate  systems  were  specified  on  connection 
bulk  data  cards  for  the  elements. 


5159  ***  USER  FATAL  MESSAGE  5159,  IN  MODULE  RETEMP  THE  SUBCASE 

NUMBER  ****  IS  GREATER  THAN  THE  NUMBER  OF  OUTPUT  TIME 
STEPS  +1,  ****. 


Subcase  numbers  after  the  first  subcase  in  Rigid  Format  13 
take  on  special  significance.  See  Section  3.14.4  of  the 
User's  Manual. 


5160  ***  USER  FATAL  MESSAGE  5160,  IN  MODULE  RETEMP  MORE  THAN  12 

RECORDS  EXIST  ON  GEOM3. 


Program  logic  error. 


5161  ***  USER  FATAL  MESSAGE  5161,  IN  MODULE  RETEMP  THE  NUMBER  OF 

WORDS  IN  AN  EQDYN  OR  SILD  RECORD  ****  IS  NOT  CORRECT. 
THE  NUMBER  OF  POINTS  IN  THE  P-SET  IS  ****. 


Program  logic  error. 
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5162  ***  USER  FATAL  MESSAGE  5162,  IN  MODULE  RETEMP  AN  SIL  NUMBER  IN 

EQDYN  IS  ****  DURING  A  BINARY  SEARCH  WE  CANNOT  FIND  IT 
IN  SILD.  THE  EXTERNAL  NUMBER  IS  ****. 

Program  logic  error. 

5163  ***  USER  WARNING  MESSAGE  5163,  IN  MODULE  STRSAVG ,  OUTPUT 

DEVICE  TYPES  FOR  ELEMENTS  ****  AND  ****  ARE  NOT  THE 
SAME; THE  OUTPUT  DEVICE  WILL  BE  THAT  FOR  THE  FIRST  ELEMENT. 

Two  elements  whose  stresses  are  being  averaged  have 
different  output  device  types,  e.g.,  print  and 
punch . 

5164  ***  USER  FATAL  MESSAGE  5164,  IN  MODULE  STRSAVG,  SCALAR  INDEX 

NUMBER  ****  CANNOT  BE  FOUND  WHEN  TRYING  TO  CONVERT  TO  THE 
GRID  POINT  NUMBER. 

Program  logic  error 


165  ***  USER  WARNING  MESSAGE  5165,  THERE  ARE  NO  TEMPERATURE- 

DEPENDENT  MATERIALS  IN  THIS  PROBLEM. 


5166  ***  USER  FATAL  MESSAGE  5166,  TEMPERATURE -DEPENDENT  MATERIAL 

TABLE  ****  DOES  NOT  HAVE  A  CORRESPONDING  TRANGE  CARD. 

See  the  remarks  for  the  TRANGE  bulk  data  card. 
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5167 


*** 


USER  FATAL  MESSAGE  5167,  IN  MODULE  **** 
EXIST. 


NOT  TRANGE  CARDS 


See  the  remarks  for  the  TRANGE  bulk  data  card. 


5168  ***  USER  FATAL  MESSAGE  5168,  IN  SUBROUTINE  DPDB,  THE 

NUMBER  OF  ELEMENT  ID-S  ON  A  ****  CARD  IS  NONPOSITIVE. 

Program  logic  error 


5169 


*** 


USER  WARNING  MESSAGE  5169,  DUPLICATE  ELEMENT  ID  **** 
ON  CARD  TYPE  **** 

TYPE  1  =  CONVEC 
TYPE  2  =  HGEN 
TYPE  3  =  RAD I AT I 
TYPE  4  =  RAD I AT 2 
TYPE  5  =  HFLUX 


The  duplicate  ID  is  ignored. 


5170  ***  USER  WARNING  MESSAGE  5170,  IN  A  THERMAL  TRANSIENT  ANALYSIS, 

THERE  ARE  NO  RADIAT1 ,  RADIAT2,  CONVEC,  HGEN  or  HFLUX  CARDS. 
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L 


5171 


*** 


USER  FATAL  MESSAGE  5171,  ELEMENT  ID  ****  WAS  SPECIFIED 
ON  A  CONVEC ,  RADIAT1,  RADIAT2,  OR  HFLUX  CARD  BUT  THE 
ELEMENT  IS  EITHER  IS3D8  or  IS3D20. 


Convection  and  radiation  specifications  must  be  made  for 
surfaces,  i.e.,  IS2D4,  IS2D8,  SURF4,  or  SURF8,  elements, 
not  solids. 


•5 172  ***  USER  FATAL  MESSAGE  5172,  ELEMENT  ID  ****  WAS  SPECIFIED 

ON  AN  HGEN  CARD,  BUT  THE  ELEMENT  TYPE  IS  EITHER  SURF 4  or 
SURF 8. 


SURF4  arid  SURF8  elements  are  dummy  surface  elements.  Heat 
generation  specifications  must  be  made  for  only  IS2D4, 
IS2D8,  IS3D8 ,  or  IS3D20  elements. 

5173  ***  USER  FATAL  MESSAGE  5173,  SUBROUTINE  ****  CANNOT  FIND  ANY 

IS2D4 ,  IS2D8 ,  IS3D8,  IS3D20,  SURF 4  or  SURF 8  ELEMENTS  ON 
EST. 


Program  logic  error,  or  a  thermal  transient  problem  does 
not  contain  any  of  these  elements. 


5175  ***  USER  FATAL  MESSAGE  5175,  IN  SUBROUTINE  DPDB  WE  CANNOT  FIND 

THE  REQUIRED  NUMBER  OF  ELEMENT  ID  MATCHES  OF  THE  ID-S  ON 
CARD  TYPE  ****  WITH  THOSE  ON  THE  SCRATCH  FILE. 


Check  that  there  exist  elements  with  ID-S  that  match  the  ID-S 
on  the****  cards.  See  message  5169  for  the  card  type 
correspondence. 
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5176  ***  USER  WARNING  MESSAGE  5176,  ELEMENT  ID  ****  REFERENCE  A  8-D 

ELEMENT  BUT  THE  ID  APPEARS  ON  A  CARD  TYPE  ****.  THIS 
ID  WILL  BE  IGNORED. 

CONVEC ,  RADI ATI,  RADIAT2,  AND  HFLUX  may  reference  only 
IS2D4,  IS2D8,  SURF4  or  SURF8  elements.  See  message  2169 
for  the  card  type  correspondence. 


5177  ***  USER  WARNING  MESSAGE  5177,  POINT  ID  ****,  COMPONENT  ID 

****  WERE  SPECIFIED  ON  A  TIC  CARD.  THE  COMPONENT  MUST 
BE  1  IF  THERE  ARE  NO  SCALAR  POINTS  IN  THE  PROBLEM  AND 
MUST  BE  0  OR  1  IF  THERE  ARE  SCALAR  POINTS.  THIS  TIC  WILL 
BE  IGNORED. 

In  thermal  transient  problems,  the  only  degree  of  freedom 
is  temperature  which  is  assumed  to  be  component  1. 

5178  ***  USER  FATAL  MESSAGE  5178,  NO  MATCH  ON  EQEXIN  FOR  POINT  ID 

****  ON  A  TIC  CARD. 

Obscure  program  error,  or  the  grid  point  on  a  TIC  card  was 
not  defined  on  a  GRID  card. 

5179  ***  USER  FATAL  MESSAGE  5179,  PIVOT  FROM  GPCT  ****  DOES  NOT 

MATCH  PIVOT  FROM  CONVPT. 


Program  logic  error. 
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5180  ***  USER  FATAL  MESSAGE  5180,  SIL  ****  CANNOT  BE  FOUND  IN 

GPCT  RECORD  WITH  PIVOT  ****. 

Program  logic  error. 

5182  ***  USER  FATAL  MESSAGE  5182  SCALAR  INDEX  NUMBER  ****  IS 

EITHER  IN  BOTH  THE  D-SET  AND  S-SET  OR  IS  IN  NEITHER. 


No  thermal  multi-point  constraints,  extra  points,  or 
omitted  coordinates  are  used  in  thermal  transient 
analysis . 


5183  ***  USER  FATAL  MESSAGE  5183,  IN  TRDTEM ,  THE  NUMBER  OF  DEGREES 

OF  FREEDOM  FROM  TRL  IS  ****  BUT  THE  NUMBER  OF  ROWS  OF 
THE  REDUCED  LOAD  VECTOR  IS  ****.  THESE  SHOULD  BE  EQUAL. 

Program  logic  error. 

5184  ***  USER  FATAL  MESSAGE  5184,  IN  SUBROUTINE  ****  NOT  ENOUGH 

MATCHES  BETWEEN  EST  AND  RAD1D. 

In  computing  loads  as  specified  on  RADIAT1  or  RADIAT2  cards, 
not  enough  matches  were  made  in  finding  the  element  ID-S. 

At  this  point,  this  is  probably  a  program  logic  error. 
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5185  ***  USER  WARNING  MESSAGE  5185,  IN  MODULE  TRDTEMP ,  NOT  ENOUGH 

TIME  REMAINS  TO  COMPLETE  THE  SOLUTIONS. 

See  the  explanation  to  message  3045. 

5186  ***  USER  FATAL  MESSAGE  5186,  IN  SUBROUTINE  SOLVTI ,  THE  ORDER 

OF  ROW  POSITIONS  IS  NOT  STRICTLY  INCREASING. 


Program  logic  error. 


5187  ***  USER  FATAL  MESSAGE  5187,  IN  SUBROUTINE  SOLVTI,  THERE  IS 

A  NULL  COLUMN  IN  THE  COEFFICIENT  MATRIX  OF  THE  THERMAL 
TRANSIENT  ANALYSIS. 


The  coefficient  matrix  is  singular. 
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THEORETICAL  MANUAL  ADDITIONS 


5.13  ISOPARAMETRIC  ELEMENTS 


5.13.1  Introduction 


If  the  shape  functions  chosen  to  describe  the  curvilinear 
coordinates  of  a  finite  element  are  identical  to  those  used  to 
prescribe  the  displacement  function  variation,  then  the  element 
is  termed  isoparametric. 

This  definition  and  the  formulation  of  the  elements  in  this 
section  are  derived  and  taken  from  the  formulation  described  in 
references  1,  2  and  3,  and  4.  The  isoparametric  elements  in 
NASTRAN  are: 

1.  IS2D4  -  Planar  quadrilateral  element  with  zero  bending  stiffness 

2.  IS2D8  -  Planar  quadriparabolic  element  with  zero  bending  stiffness 

3.  IS3D8  -  Solid  linear  element 

4.  IS3D20  -  Solid  parabolic  element 

The  isoparametric  family  of  elements  uses  a  mapping  technique 
to  avoid  performing  the  required  numerical  integrations  on  arbitrarily 
shaped  elements.  The  displacement  function  is  assumed  in  relation 
to  an  element  with  a  simple  shape,  e . g .> rectangle  in  two  dimensions, 
rectangular  parallellepiped  in  three  dimensions.  However,  the  shape 
of  the  element  is  permitted  to  distort  and  take  another  shape,  dictated 
by  the  same  functions  that  describe  the  displacement  patterns.  Figure 
1  shows  the  four  available  elements  and  the  simple,  or  parent,  element 
is  used. 

The  grid  points  of  each  parent  element  are  defined  in  a  local  (£, 
h,C)  coordinate  system.  The  grid  points  at  the  corners  have  values 
of  +_  1 .  The  grid  points  of  the  curvilinear  element  are  defined  in 
the  basic  (x,  y,  z)  coordinate  system.  Each  parent  element  has  an 
assumed  displacement  function,  which  is  also  used  to  relate  the  (£',n,0 
coordinates  to  the  (x,  y,  z)  coordinates. 
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Curvilinear 
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Figure  1.  Isoparametric  Elements  and  Their 
Associated  Parent  Elements 
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5.13.2  Structural  Elements 

5.13.2.1  Planar  Elements 

Consider  the  planar  elements  in  Figures  2a  and  2b  below. 
(The  order  of  the  grid  point  numbering  is  required.) 


Figure  2b.  IS2D8  Element 
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These  elements  have  a  thermal,  as  well  as  structural,  capability 
and  may  have  anisotropic  thermal  and  structural  temperature-dependent 
material  properties. 

For  the  thermal  and  structural  elements,  lumped  or  consistent 
capacitance  and  mass  matrices,  respectively,  may  be  computed. 
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Let  u  and  v  be  the  components  of  displacement  in  the  element 
coordinate  system,  parallel  to  the  x-  and  y-directions  of  the  element 
coordinate  system,  respectively.  The  in-plane  displacements  of  all 
the  grid  points  of  the  element  are  represented  by  the  vector 


/ 


<V‘  1 


a) 


111 

where  n=  4  or  8  depending  on  the  element.  If  [K  ]  is  the  associated 

Vv 

stiffness  matrix,  then 

[Kee](ue)=(feJ  (2) 

where  (f  }  is  the  vector  of  in-plane  forces, 
e 

[K  1  is  derived  as  follows: 

Assume  that  the  displacement  function  is  given  by 

{u}= [N]  {ue}  (3) 

where  {u}  is  the  vector  of  in-plane  displacements  and  [N]  is  tjie  shape 

matrix.  The  shape  matrix  [Nj  may  be  given  as  [N]  =  [IN^, 

where  I  is  the  2x2  identity  matrix.  The  are  defined  as  follows: 

IS2D4 

Ni  =  j(l+C5i)(l+nni)  (4) 


IS2D8 

(i)  Comer  grid  points 

(ii)  5**1,  n=0 

(iii)  5=0,  tf*! 


Ni  =  jU+ey  (1  +  nn.)  (se.+wyl) 

N.=  i(l+55.)(i  -n2) 

Nk  =  j(l-52) (l+nnk) 

(5) 


where  i,  j,  k  take  on  values  of  the  grid  point  numbers. 
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The  membrane  strains  are  given  by: 


where  the  x  and  y  subscripts  refer  to  partial  derivatives  of 

with  respect  to  the  x-and  y- directions  of  the  element  coordinate 
system.  Note  that  these  partial  derivatives  cannot  be  directly  ob¬ 
tained  since  the  N^  functions  are  functions  of  the  (£,n  )  coordinate 
system.  Therefore,  we  must  define  the  Jacobean  fj]  as  follows: 
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where 


t  is  the  element  thickness 

G  is  the  material  properties  matrix 
6 

det  [J]  is  the  determinant  of  the  Jacobean. 

(It  can  be  shown  that  det  [J]  d£dn=dxdy) . 

In  NASTRAN  the  [Gg]  matrix  may  be  anisotropic  so  that  the  only 
restriction  is  that  it  be  symmetric.  In  the  case  of  isotropy. 


[Ge]: 


1-v 

i>E 

l  V 


vE 

lV 

E 

lV 


(13) 


where  the  shear  modulus  G  is  G  =  E/2(l+vj. 

If  anisotropic  materials  are  used,  the  user  specifies  them  with 
respect  to  a  particular  orientation,  which  does  not  necessarily 
correspond  to  the  principal  axes.  The  angle  6  (see  figures  2a  and 
2b)  is  an  input  parameter.  The  properties  matrix  in  the  material 
coordinate  system  is  transformed  into  the  element  coordinate  system 
by 


[ge]  =  [u]t[gj  [U] 


where 


[U]  = 


„  2  n 

.  2 

COS  0 

sm  0 

cos0sin0 

.2 

2 

sin  e 

cos  0 

-cos0sin0 

-2cos0sin0 

2cos0sin0 

2  .  2 
cos  0-sm  0 

(14) 
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is  the  transformation  matrix  for  the  rotation  of  strain  components. 


The  integration  specified  in  equation  (12)  is  performed  numerically 
using  Gaussian  quadrature.  The  stiffness  matrix  is  approximated  by 


m=l 


(15) 


where 


A  is  the  coefficient  used  in  the  quadrature  and  depends 
on  the  number  of  quadrature  points 


p  is  the  number  of  quadrature  points 
f  =  [B]T[Ge][B]det[J] 


The  number  of  quadrature  points  p  used  for  the  IS2D4  element  is  2  for 
IS2D8,  2  or  3.  The  number  p  may  vary.  As  p  increases,  theoretically, 
the  solution  improves.  However,  running  times  increase  rapidly  as  p 
increases.  The  values  of  p  used  for  these  elements  appear  to  be  very 
satisfactory. 

Finally,  the  stiffness  matrix  is  transformed  from  the  local 
element  coordinate  system  to  the  global  coordinate  system  of  the  grid 
points.  If  this  transformation  is  given  by 

{ug}  =  [T] (Ug) ,  (16) 


then 

OCGG>  »  mTIKee][T]  (17) 

The  consistent  mass  matrix  for  these  planar  elements  is  given  by 
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[Mee]  =  t;-l  ;-l  det  [J]d?dr<  (18) 

where  p  is  the  mass  density.  The  appropriate  transformation 
gives  [Mgg] . 

The  lumped  mass  matrix  is  formed  by  distributing  the  mass 
proportionately  among  the  grid  points  of  the  element.  This  matrix 
is  diagonal.  The  total  mass  for  an  element  is  PV,  where  V  is  the 
element  volume,  so  that  the  mass  applied  to  grid  point  i  is 

/nJ 

~  Pt/_1  f_1  det[j]d£dp.  (19) 

E  / N2 
j-1  J 

A  load  vector  is  produced  by  the  thermal  expansion  of  an  element. 
The  thermal  strain  vector  is 


{et}= 


where  (am>  is  the  thermal  expansion  coefficient  vector  and  is 
the  temperature  at  element  grid  point  i,  and  TQ  is  the  reference 
temperature. 


Then, 


So,  the  load  is 


<v-  -  t<y<WTo>  • 


{Pe>=  /‘j  [B]Tdet[J]dEdri]  [Ge]{om}(ti-T0)  (22) 


Transformation  to  global  coordinates 


{PG)=[T]1{Pe> 


completes  the  calculation. 
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After  the  grid  point  displacements  have  been  computed,  element 


stresses  are  computed  as  follows: 

(ueMT]  {ug} 

t — \ 

(M 

V - ' 

{e->=[B]{ue} 

(25) 

{o}=[GM]{e-et} 

(26) 

or 

-v 

(27) 

By  evaluating  [B]  at  the  gauss  quadrature  points,  stresses  are 
computed  at  these  gauss  points.  Grid  point  stresses  are  computed 
by  extrapolating  from  the  gauss  points.  This  method  appears  to  be 
superior  to  computing  grid  point  stresses  directly.  Also,  the 
number  of  gauss  quadrature  points  used  in  computing  the  mass 
matrix  is  one  more  than  that  used  in  computing  the  stiffness 
matrix. 
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5.13.2.2  Solid  Elements 

The  two  solid  elements  in  NASTRAN  are  the  IS3D8  and  IS3D20  ele 
ments,  which  are  connected  to  8  and  20  grid  points,  respectively. 
The  sides  of  the  IS3D8  element  are  linear;  the  sides  of  the  IS3D20 
element  are  parabolic.  An  IS3D8  element  is  shown  in  Figure  3. 


z 


Many  details  in  the  derivations  below  are  very  similar  to  those 
for  the  planar  elements  and  will  not  be  repeated. 


If  anisotropic  material  properties  are  to  -be  specified  for  a 
solid  element,  the  user  supplies  a  coordinate  system  number  which 
refers  to  a  rectangular  material  coordinate  system.  The  element 
geometry  is  then  transformed  into  this  material  coordinate 
system,  in  which  all  calculations  are  made  until  the  final 
conversion  to  the  global  coordinate  system. 

Let  u,  v,  w  be  the  components  of  displacement  in  the  material 
coordinate  system,  parallel  to  the  x-,  y-,  z-axes  of  that  system, 
respectively.  Then  the  displacements  of  the  grid  points  of  the  ele 
ment  are  represented  by  the  vector 
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<ueH 
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where  n  is  8  or  20,  depending  on  the  type  of  element.  The  shape 

functions  N.  are  as  follows: 

1 


IS3D8 

N.=|(l+£q)(l+nni)(l+^i)  (2) 

IS3D20 


(i) 

Corner  grid  points 

(l+nni)  (l+C^)  (5?.+nn 

(ii) 

xTH 

» — i. 

II 

O 

N  4(l-52)(l+nn.)(l+ccJ 

J  4  3  3  (3) 

(iii) 

v° 

Nk4(l+5Ck)(l“n2)(l+«k) 

(iv) 

V° 

The  shape  matrix  [N]  may  be  represented  as 


[N]=[IN1,  IN2,...INn] 


(4) 


where  I  is  the  3X3  identity  matrix.  Then  the  displacement  function 
is  given  by 


(u)  =[N]{ue> 


The  strains  are  given  by 


(e}=. 
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(5) 


where 
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The  Jacobean  [J]  becomes 
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The  stiffness  matrix  [K  ]  then  becomes 

G  G 

[Kee]=/il  f\  [B]T[Gm][B]det[J]dCdnd?  (8) 

where  [Gm]  is  the  material  properties  matrix  which,  in  the  case  of 
isotropy,  is: 


1-v  1-v 
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The  approximation  for  the  integration  is 


[K  ]* 

L  eeJ  Jt=l  m 


II  jlAA,Vf«l’V'n)l 


(10) 


For  IS3D8  elements  p=2,  for  IS3D20  elements  p=3.  Finally,  if 
[T]  is  a  transformation  matrix  from  the  material  coordinate  system 
to  the  global  coordinate  system,  then 

[KGGMT]T[Kee][T]  (11) 


The  consistent  mass  matrix  is  given  by 

[Mee]=P/-l  ;-l  ;-l  [N]T[N]det[J]d?drld^ 

The  lumped  mass  matrix  is 
1  1  1 

[M  ]=p/  /  .  f  ,  det[J]d?dndc 

ee  -1  -l  -l 

and  the  mass  associated  with  grid  point  i  is 

/N?  Ill 

— - 1  det[J]d5dndc 

j=l  1 

The  load  vector  due  to  thermal  expansion  is 

ft i  (14) 

For  completeness,  the  element  stresses  are,  again, 
(o]=[Gm][B][T](ug>-[GJ{c.m}(trT0) 

Grid  point  stresses  are  computed  by  first  computing  stresses  at 
gauss  quadrature  points  and  then  extrapolating  to  the  grid  points. 

Also,  the  number  of  gauss  quadrature  points  used  in  computing 
the  mass  matrix  is  one  more  than  that  used  in  computing  the  stiffness 
matrix. 


(12) 


(13) 
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5.13.3  Thermal  Elements 

The  differential  equation  governing  the  behavior  of  the  tem¬ 
perature  u  is 


3,,  3u..  3  ,,  3u.  3..  3u.  r3u  n 

ta(kx  =0 


(0 


where 


G 

C 


are  the  anisotropic  conductivity  coefficients 
is  the  rate  of  heat  generation  or  external 
heat  input,  and 
is  the  capacitance 


In  addition  to  specified  temperatures  on  the  boundaries,  another 
boundary  condition  might  take  the  form 


(2) 


where 


q 

H 

c 


are  the  direction  cosines  of  the  outward  nor¬ 
mal  to  the  boundary  surface 
is  the  heat  flux  per  unit  of  surface 
is  the  film  coefficient  for  convection  ,  and 
is  the  ambient  temperature  for  convection 


If  k^,  ky,  k^  are  equal  and  q  and  a  are  zero,  then  equation  (2) 
reduces  to  a  condition  applicable  to  non-conducting  boundaries 


3u 

3n 


=0 


(3> 


By  applying  Euler's  theorem  in  the  calculus  of  variations,  equation 
(1)  will  be  satisfied,  subject  to  boundary  condition  (2),  if  and 
only  if  the  functional 

[I{kx O 2+ vfy’ 2*kz  O 2) '  (g'c!t)u1 m 

is  minimized  and  if  the  boundary  conditions  (2)  are  satisfied. 
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Zienkiewicz1 states  that  boundary  conditions  (2)  may  be  taken  into 
account  by  modifying  equation  (4)  as  follows: 

t?  (kx02tky'f7’2*kz(l7)2-(G-C!r)u]d)'d>'dz 

+/D (q-H  u  )udS+/D  u2dS  (4 1 ) 

c  a  B  2  c  k  J 

where  the  last  two  integrals  are  surface  integrals  taken  over  the 
boundary  on  which  conditions  (2)  are  applied.  (We  will  not  be 
too  careful  about  the  placement  of  the  thermal  coefficients  H  , 
a,  k^,  . ..,  etc.  since  these  parameters  are  considered  to  be  con¬ 
stant  within  an  element.) 

Now  assume,  as  in  previous  sections,  that,  for  any  element, 

u=[N]{ue>  (5) 


where  now. 


[N]=[N1,  N2,  ...,  Nn], 


,  and 


n  is  the  number  of  grid  points  associated  with  the  element 
then  /  \ 


u=[N1>  N2,  ... .N^] 


u] 

u„ 


> 


/ 


(6) 


We  now  want  to  minimize  the  functional  X  with  respect  to  all  u^’s 

in  the  problem.  We  can  do  this  by  evaluating  the  contributions 

3X 

to  each  differential,  g^j~,  from  each  element,  adding,  and  equating 
to  zero.  If  the  contribution  to  X  from  an  element  is  Xg,  then, 
differentiating  equation  (4’) 
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(7) 


e  -fifty  9u  9  (9  Hk  —  —  r— i+k  —  —  (— ) 

9ir  ~fff[  x  97  30TC3F  ky  a-  ^  Kz  *" 

l  i 


3y  3ui  l3yJ  z  3z  3^  l3zJ 


-(g-cIt)  It1  dxdydz 

i 


t/Jq-Hu)  dS+/nH_u  dS 


c  a'  3u. 

l 


Be  3u. 

l 


Now,  substituting  equation  (6)  into  equation  (7),  we  obtain 


3X 


~  =///{k  N.  fN,  N  ]+k  N.  [N.  ]+k  N.  [N, 

3u.  x  ixL  lx  nxJ  y  iyL  ly*  ’  nyJ  z  izL  lz 


dxdydz{ue) 

3u 


-///CN  dxdydz+///CNi[N1, . .  .Njdxdydz  ~ 

+/D( q-H  u  )N. dS+/nH  N.[N.,  N  ]dS{u  } 

B  n  c  a  i  Beil  nJ  e 

If  all  the  contributions  from  an  element  are  denoted  by 


(3*  } 

1 3u  se 


I 


3X 
_ £ 

3u, 


3X 

£ 

3u 


n  / 

then  the  notation  can  be  reformulated  to 


,  N  ]} 
’  nzJ 


(8) 


(9) 


{—  }  =[K  ]{u  }+[B  ] 

3u  e  L  eeJ  e  L  eeJ 


3u 

{ — -  } 
1  3t 


<Pe> 


(10) 


where 

[K  ]  is  the  element  conduc  tivity  (stiffness)  matrix 

[B  ]  is  the  element  capacitance  (damping)  matrix,  and 
6  6 

{Pg  }  is  the  vector  of  element  loads. 


Equating  equation  (10)  to  0  and  performing  some  manipulation  yields, 
for  the  entire  problem, 
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[B] (u)+[K] {u}  =-{P) 

ni) 

Then, 

,=zll/(kN.  N.  +k  N.  N.  +k  N.  N.  )dxdydz+E/DH  N.N.dS 

J  vx  ix  jx  y  ly  jy  iz  jzJ  1  B  c  l  3 

(12) 

B.  .  =£///CN . N .dxdydz 

13  i  J 

(13) 

and 

P.=-E//JCN  . dxdydz+E/D (q-H  u  )N.dS 
l  l  1  Bvn  c  a  i 

(14) 

where  the  summation  covers  the  contributions  from  each  element. 
Finally,  the  matrices  and  vectors  are  assembled  in  exactly  the 
same  way  as  in  the  structural  problem. 
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19.  THERMAL  TRANSIENT  -  STRUCTURAL  STATICS  COMBINATION 


Thermal  transient  -  Structural  Statics  Combination,  Rigid  Format  14, 
is  a  combination  of  the  analyses  contained  in  Direct  Transient  Analysis, 
Rigid  Format  9,  with  modifications,  and  Static  Analysis,  Rigid  Format  1. 

In  the  transient  analysis  portion,  the  temperature  history  of  the  structure 
is  computed,  taking  account  of  such  thermal  conditions  as  radiation, 
convection,  flux,  and  heat  generation.  Then,  a  series  of  static  analyses 
will  be  performed  using  ,as  part  of  the  static  loading  the  equivalent 
loads  due  to  the  temperature  distribution  at  times  selected  by  the  user. 

The  finite  elements  allowed  to  be  used  in  Rigid  Format  14  are  the 
isoparametric  elements  IS3D4,  IS3D8,  IS3D8,  and  IS3D20,  the  surface 
elements  SURF1,  SURF4,  and  SURF8,  and  the  BAR  element.  These  elements 
may  have  anisotropic,  temperature-dependent  material  properties  in  both 
the  thermal  and  structural  analyses.  Therefore,  in  the  thermal  transient 
analysis,  element  conductence  and  capacitance  matrices  may  be  recomputed. 
Recomputation  of  element  matrices  will  also  be  performed  if  convection  is 
specified. 

A  very  simplified  flow  diagram  of  Rigid  Format  14  is  given  in 
Figure  1. 
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19.1  THERMAL  TRANSIENT  ANALYSIS 

In  thermal  transient  analysis,  the  governing  differential 
equation  is  a  first  order  equation  rather  than  the  second  order 
equation  governing  the  structural  problem.  Before  the  integration 
algorithm  is  described,  the  thermal  transient  loads  will  be  described. 

19.1.1  Thermal  Transient  Loads 

The  following  types  of  loading  are  available  in  a  thermal 
transient  analysis: 

1.  Transient  loads  of  the  types  discussed  in  Section  11.1 

2.  Nonlinear  loads  of  the  types  discussed  in  Section  11.2 

3.  Convection  from  a  surface  to  its  surroundings 

4.  Radiation  from  a  surface  to  its  surroundings 

5.  Radiation  from  one  surface  to  another  surface 

6.  Internal  heat  generation 

7.  Heat  flux 

Of  the  new  types  of  loading,  convection,  heat  generation, 
and  heat  flux  were  discussed  in  Section  5.13.3.  An  expansion 
of  those  remarks,  as  well  as  a  discussion  of  the  radiation  loads, 
will  be  presented  here.  (Equations  (12)  (13),  and  (14)  in  Section 
5.13.3  include  convection,  heat  generation,  and  heat  flux.) 

Convection  loads  are  given  by 

Q=HcA(T-TJ  (1) 

where 

Q  is  the  heat  flow  per  unit  time  due  to  convection  exchange 

Hc  is  the  film  coefficient,  which  may  be  time  dependent 
A  is  the  area  of  the  surface 
T  is  the  temperature  of  the  surface 
T^  is  the  ambient  temperature 

Contributions  of 

H  //N.N.dS  (2) 

cJ'  i  j  v  1 

are  made  to  the  (i,j)th  entry  of  the  conduction  matrix  for  the 
element.  The  contribution  to  the  transient  load  vector  from  each 
grid  point  defining  the  surface  is 

HJJVN-dS  (3) 
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where  and  fL  make  up  the  element  shape  matrix  [N] .  (See  Section 
5.13.2.1.)  (The  integrals  in  equations  (2)  and  (3)  are  surface 
integrals,  as  are  ffds  integrals  in  equations  (8),  (10),  and  (12).) 
Internal  heat  generation  loads  are  given  by 

Q=GV  (4) 

where 

Q  is  the  load  to  be  applied 

G  is  the  heat  generation  per  unit  time,  which  may  be  tempera¬ 
ture  dependent,  and 
V  is  the  element  volume 

For  solid  elements,  the  contribution  from  each  grid  point  of  the 
element  is 


GCTi)///Nidxdydz  (5) 

while  for  planar  elements,  the  contribution  is 

G(Ti>h//N.dxdy  (6) 

where  h  is  the  element  thickness  and  is  the  temperature  of  grid  point 
Loads  due  to  heat  flux  are  given  by 

Q=qA  (7) 

where 

Q  is  the  load  to  be  applied 

q  is  the  heat  flux  coefficient,  which  may  be  time  dependent, and 
A  is  the  area 

The  contribution  to  the  load  vector  from  each  grid  point  defining 
the  surface  is 


//qCt)N.dS 


(8) 
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The  load  due  to  radiation  from  a  surface  to  its  surroundings 


Q=HA  (T4 -T^  4  )  (9) 

where 

Q  is  the  net  rate  of  radiation  heat  exchange  between  a 
specified  surface  and  its  surroundings 
H  is  the  radiation  factor,  which  may  be  temperature  dependent 
A  is  the  area  of  the  surface 
T  is  the  temperature  of  the  surface,  and 
Too  is  the  ambient  temperature,  which  may  be  time  dependent. 

The  contribution  to  the  load  vector  from  a  grid  point  defining 
the  surface  is 

~H(T  ) (T4-T4)//n  dS 

1  1  (10) 

where  is  the  temperature  of  grid  point  i. 

For  radiation  between  two  surfaces, 

q=ha1a2(t14-t24)  (11) 


where 

Q  is  the  net  rate  of  radiation  heat  exchange  between  surfaces 
Sj  and  S2 

H  is  the  radiation  factor,  which  may  be  temperature  dependent 
A^,A2  are  the  areas  of  and  S 2,  respectively,  and 
T1,T2  are  the  temperatures  of  and  S2>  respectively. 


The  contribution  to  the  load  vector  from  a  grid  point  defining  surface 
is 


■H(T.)(T14-T24)A2//N.dS 


(12) 


and,  from  surface  S2,  the  contribution  is 


(13) 
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(In  equations  (2),  (3),  (8),  (10),  (12)  and  (13),  the  surface 
integral  S=//dS  is  computed  as  follows. 

S=//^G-F2d?dn  (14) 


„_,3x.2  r3y,  2  ,3z.2 

E=(3e°  +  (3c° 


3x  3x  3^_  3^  3z  3z 
35  3n  *  !{  3n  35  3n 


li,2.  f|X!2  +  fe2 


and 


(Notation  is  the  same  as  in  Section  5.13.S) 


(16) 

(17) 

(18) 


19.1.2  Integration  Algorithm 

The  system  of  differential  equations  to  be  solved  is 

[Bg|  +  K](u}={P}  (1) 

Although  the  coefficient  matrices  [B]  and  [K]  may  be  functions  of 

time  t  and  of  the  solution  vector  (u),  it  is  assumed  that  they 

vary  slowly  so  that  they  can  be  considered  constant  during  any 

one  time  step.  The  load  vector  {  P)  may  be  the  sum  of  a  time  dependent 

load  {TL}  and  a  non-linear  load  {P.J. 

T  N 

The  solution  algorithm  begins  with  initial  conditions  pre¬ 
scribed  at  time  ,  and  computes  the  solution  at  prescribed  times 
t^  i=l ,  2,  3,.  .  .. 

If  at  time  t^,  i=0,  1,  2,  .  .  .,  the  solution  {u^}  is  prescribed 
or  has  been  previously  computed,  the  algorithm  proceeds  to  compute 
the  solution  {u^+j}  at  t^ie  following  time  t^+^.  For  the  time 
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increment  At^=t^+j-t^,  t^ie  solution  of  Equation  (1)  is 

{ui+l>  =  /ttl+1  {P)dt  +  exp(-At[B_1K]){u.}  .  (2) 

i 

This  is  the  value  of  the  solution  to  be  approximated. 


The  Approximation  of  the  Solution  {u} 

Since  the  load  vector  is  the  sum  of  a  time  dependent  load  {PT> 
that  can  be  computed  at  t^  and  t^+j,  and  a  non-linear  load  that  is 
computed  only  at  t^,  the  integral  in  Equation  (2)  is  approximated 
by 


At{  i(P_  +PT  )+PN  }  . 
i  i+1  i 


(3) 


The  matrix  in  exponential  form  from  Equation  (2)  in  a  series 
expansion, 

[I]  -At[B-1K]  +  [  B- 1 K]  2  -  [B_1K]3+...,  (4) 

2 

is  equal,  up  to  terms  of  order  (At)  ,  to  the  expression 


[B+  ^  K]  1 [B-  K]  =  [1+  B-1K]  1 [I-  ^|b_1K] 

=  [I]  -  At[B-1K]  +  [B_1K]  2-  [B~1K]  3+  -  •  •  . 


(5) 


Thus  an  approximate  solution  in  the  form  of  the  solution  of 
a  system  of  linear  algebraic  equations  is  given  by 


[B+  ~  K]  {u.  +  1}  =  [B-  K]  {u.}  +  [B+  ^  K]{i(PT  +PT  )+PN  }  At 

i  i+1  Ni 

(6) 

This  linear  system  is  solved  by  a  symmetric  decomposition  and  forward 
backward  substitution. 
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